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I. 


SUMMARY 


Through  experimental  and  theoretical  studies,  the  Penn  nonlinear  optics  research 
program,  principally  funded  by  DARPA  and  AFOSR,  has  firmly  established  that  the  largest 
second  ("^>>31®!  1^2)  third  nonlinear  optical 

properties  occur  in  organic  and  polymer  systems.  This  major  development  together  with 
progress  in  optical  bistability,  phase  conjugate  wave  generation,  and  neural  networics  has 
stimulated  considerable  research  and  development  activities  in  centers  throughout  the  world 
with  many  potential  applications  being  pursued  in  optical  signal  processing,  integrated 
optics,  image  reconstmction,  sensor  protection,  data  storage,  and  optical  memory  and  logic 
technologies.  Most  notably,  many  industrial  l^x>ratories  have  successfully  developed 
electrooptic  polymers  and  devices  based  on  the  experimental  and  theoretical  results 
achieved  by  the  current  research  program. 

The  major  purpose  of  the  research  program  is  to  further  rapidly  accelerate  and 
expand  these  advances  in  nonlinear  optics  research.  The  principal  approach  and  objectives 
during  the  last  grant  period  were  comprised  of  the  following  original  statements, 
o  Investigate  the  physical  mictoscopic  origin  and  mechanism  of  second  and  third 
order  nonlinear  optical  responses  in  conjugated  polymers, 
o  Develop  new  theoretical  concepts  relating  to  the  many-body  nature  of  highly 
charge-correlated  rc-electron  states. 

o  Calculate  nonlinear  optical  responses  of  new  71-electron  topological  structures, 

o  Determine  frequency  dependent  third  order  and  112  respones  at  fast  time  scales, 

o  Demonstrate  optical  bistability  in  nonlinear  optical  conjugated  polymers. 

The  remainder  of  the  report  is  organized  in  two  major  sections.  Section  II  contains 
reviews  of  Selected  Research  Highlights  from  the  research  program  during  the  past  two 
year  funding  period.  Section  HI  provides  an  extended  description  of  major  advances 
illustrative  of  the  results  obtained  by  the  funded  research  program.  This  section  contains 


the  subjects  (I)  Nonlinear  Optical  Processes  in  Lower  Dimensional  Conjugated  Structures, 
and  (n)  Frequency  Dependence  of  the  Optical  Kerr  Effect  and  Third  Order  Electronic 
Nonlinear  Optical  Processes  of  Organic  systems. 
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II. 


SELECTED  RESEARCH  HIGHLIGHTS 


I.  ELECTRON  DENSITY  CORRELATION  DIAGRAMS.  MICROSCOPIC 
ORIGIN  AND  MECHANISM  OF  NONLINEAR  OPTICAL  PROCESSES  IN 
CONJUGATED  LINEAR  CHAINS 

The  fundamental  problem  of  understanding  the  origin  of  large,  nonresonant  second  order 
macroscopic  nonlinear  optical  susceptibilities  observed  for  conjugated  ;c-electron  organic  and 
polymer  solids  is  currently  viewed  as  reducible  to  theoretical  and  experimental  studies  of  the 
corresponding  microscopic  susceptibilities  of  single  molecular,  or  polymer  chain,  units  making  up 
the  nonlinear  optical  mediiun. 

Linear  polyenes  and  polyenynes  represent  the  finite  chain  limit  to  the  corresponding  quasi- 
one-dimensional  conjugated  polymer  chain,  and  phase  matched  second  harmonic  generation,  for 
example,  has  been  demonstrated  earlier  by  us  for  noncentrosymmetric  diacetylene  conjugated 
polymer  structures.  In  this  regard,  for  the  fundamentally  important  case  of  a  noncentrosymmetric 
linear  polyenyne  NMDVDA  ftrans-trans  l-(4  methylphenyl)-8-(4'-nitrophenyl)-l,  7-diene-3,5- 
octadiyne], 


we  reported  last  year  the  first  results  of  DC  induced  second  harmonic  generation  (DC  SHG) 
measurements  of  ^|j|^('2o);o),a))  (Figure  2.1.1). 

Electron  conelations  play  an  important  role  in  the  nonlinear  optical  responses  of  conjugated 
organic  structures,  and  their  description  of  microscopic  second  and  third  order  nonlinear  optical 
responses  differs  markedly  from  independent  particle  models.  As  a  result,  large  charge 
asymmetries  occur  in  the  virtual  excitation  processes,  and  we  have  long  sought  a  quantitative 
diagrammatic  description  of  these  features  of  the  fundamental  mechanism.  This  year  we 
successfully  developed  a  diagrammatic  theoretical  formalism  that  liitks  the  results  of  our  many- 
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election  calculation  of  Pjj|^(-2a):ti),a>)  by  SCF-MO-SDCI  methods  with  experimental  results 
(Figure  2.1.1). 

Election  density  correladon  diagrams  based  on  the  electron  state  functions  provide  a 
compact  graphical  representation  of  important  microscopic  features  for  second  order  nonlinear 
optical  processes.  The  transition  moment  jXpj,.  is  expressed  in  terms  of  the  transition  density 
matrix  Ppp.(r,r’)  by  ■  *0  /  dipole  moment  difference  Ap,!, 

by  the  difference  density  function  Pn'Pg  between  the  excited  and  ground  state  functions 
<Ap.p>  ■  -e  J  r  (pp-Pg)  dr  where  p^  is  the  first  order  reduced  density  matrix. 

Election  density  correlation  diagrams  Pp~Pg  and  Pp,^.  for  second  order  processes 
in  NMDVDA  involving  the  dominant  contributing  states  I,  III  and  FV  are  plotted  as 
examples  in  Figure  2.1.2.  The  solid  and  dashed  lines  correspond  to  negative  and  positive 
values,  respectively.  In  each  case,  upon  virtual  excitation  in  the  second  order  process, 
large  asymmetric  redistribution  of  electrcm  density  occurs  along  the  conjugated  chain, 
especially  in  the  x  direction,  resulting  in  large  dipole  moment  differences  App  (Figures 
2.1.2a  and  2.1.2b)  and  transition  moments  p^p.  (Figure  2.1.2c)  that  determine  P^.  Thus, 
for  example,  Piv'Pg  of  Figure  2.1.2b  corresponding  to  Ap*y  of  33.5D  exhibits  a  large 

correlation  distance  between  positive  and  negative  regions  along  the  x-direction  in 
comparison  to  Pfpg  of  Figure  2.1.2a  with  its  associated  value  of  Apj*  of  18.7D.  The 
same  charge  correlated  features  are  exhibited  by  Puy  of  Figure  2.1.2c  where  the  P|* jy  of 
12.8D  is  relatively  large.  Thus,  our  electron  density  correlation  diagrams  provide  a 
quantitative  diagrammatic  foimalism  for  nonlinear  optical  processes.  Moreover,  the 
formalism  is  general  and  can  represent  the  microscopic  features  of  second  or  third  order 
processes  for  any  conjugated  structure. 


n.  DISPERSION  OF  RESONANT  SECOND  HARMONIC  GENERATION  IN 

ORGANIC  NONLINEAR  OPTICAL  ULTRATHIN  FILMS. 


Organic  and  polymer  structures  possess  unusually  large,  ultrafast  second  and  third 
order  macroscopic  nonlinear  optical  properties  that  are  fundamentally  important  to  the  field 
of  nonlinear  optics  and  provide  previously  unachievable  design  features  for  device 
implementation  in  optical  technologies.  The  macroscopic  nonlinear  optical  properties  of 
multilayer  thin  films  can  be  formed  and  tailored  as  nonlinear  optical  ultrathin  films  by 
Langmuir-Blodgett  techniques.  Earlier,  using  novel  x-ray  scattering  and  data  analysis 
methods,  we  had  reported  the  profile  structures  of  such  films  which  possess  mosaic 
spreads  much  less  than  semiconductor  single  crystals  and  allow  sufficient  accuracy  and 
spatial  revolution  to  accurately(±l  A)  position  molecular  bonding  sequences  in  a  single 
bilayer. 

In  this  report,  we  describe  the  first  experimental  and  theoretical  results  of  resonant 
second  harmonic  generation  of  a  single  asymmetric  organic  monolayer  of  a  cyclic  polyene 
DCNQI  (Figure  2.2.1).  Figure  2.2.2  shows  the  frequency  dependent  microscopic  second 
order  nonlinear  optical  susceptibility  Pj||^(-2a);cD,0))  of  a  DCNQI  monolayer.  Included  are 

off-resonance  results  obtained  from  separate  DCSHG  liquid  solution  measurements.  With 
increased  frequency,  the  values  starting  at  approximately  -30x10'  esu  increase 

smoothly  to  a  resonant  two-photon  maximum  of  -260x10'  esu  near  800  nm 
corresponding  to  the  singlet-singlet  excitation  to  the  first  p-electron  excited  state  of  DCNQI. 

The  theoretical  solid  curve  in  Figure  2.2.2  was  calculated  by  SCF-MO-SDCI 
methods  in  the  rigid  lattice  CNDO/s  approximation  that  included  both  singly  (SCI)  and 
doubly  excited  configurations  (DCI)  to  account  for  electron  correlations.  Comparison  of 
theory  and  experiment  is  quite  satisfactory. 

Density  matrices  of  the  state  functions  provide  a  compact  graphical  representation  ot 
important  microscopic  features  for  second  order  nonlinear  optical  processes.  Contour 
diagrams  for  the  states  principally  responsible  for  are  plotted  in  Figure  2.2.3  where  the 


solid  and  dashed  lines  correspond  to  negative  and  positive  values,  respectively.  In  each 
case,  upon  virtual  excitation  in  the  second  order  process,  large  asymmetric  redistribution  of 
electron  density  occurs  with  large  correlation  distances  across  the  conjugated  structure, 
especially  in  the  large  dipole  x  direction.  This  correlated  Jt-electron  redistribution  results  in 
large  dipole  moment  differences  and  transition  moments  that  determine 
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Figure  2.2.3 


III.  DESIGN  OF  NOVEL  STRUCTURES  WITH  UNUSUALLY  LARGE  SECOND 
ORDER  RESPONSES 

It  is  now  widely  accepted  that  conjugated  organic  polymer  materials  possess  the 
largest  electronic  nonlinear  optical  responses.  Use  of  organic  materials  in  optical  device 
implementation  has  already  been  under  intense  development  in  several  industrial 
laboratories  and  considerable  success  has  been  achieved,  particularly  in  the  production  of 
electrooptic  guest-host  polymer  thin  film  devices.  To  further  improve  the  performance  of 
these  devices,  it  has  now  become  more  and  more  important  that  specific  molecular 
structures  with  larger  second  order  susceptibilities  3jj|^(-(03;co^  ,0)2)  be  developed.  We 
have  long  held  the  view  that  in  order  to  achieve  advances  in  the  design  of  Ji-clectron 
structures  with  large  values  of  3jj|((*<D3;(D^  ,(02),  it  is  essential  to  have  a  complete 
microscopic  understanding  of  second  order  excitation  processes  in  these  many-electron 
systems.  To  this  end  we  have  successfully  developed  a  many-cleciron  formalism  and 
method  for  the  calculation  of  microscopic  second  order  nonlinear  optical  responses.  The 
results  of  this  method  provides  direct  insight  and  guidelines  toward  the  design  of  new 
molecular  systems. 

The  molecular  dye  MFC  [N-methyl[4(lH)-pyridinylidenc  cthylidene]2,5- 
cyclohexadien-l-one]  has  been  reported  by  a  French  research  group 

—'■1 

to  have  a  measured  Pj,(-2(o;(0,a))  value  of  more  than  1000x10’^^  esu  far  off  resonance. 

However  we  have  studied  MFC  theoretically,  and  found  by  SCF-MO-SDCI  methods  that 
the  calculated  value  was  much  smaller  by  an  order  of  magnitude.  On  the  basis  of  our 
calculated  results,  we  initiated  detailed  experimental  studies  of  3,^  at  our  laboratory  together 


with  collaborative  studies  at  AT&T  labcnatories  that  acted  as  blind  tests.  The  results  from 
both  measurements  confirmed  not  only  the  calculated  value  of  ^ut  also  the  sign  and 

frequency  dependence  (Figure  2.3.1).  The  SCF-MO-SDCI  method  indeed  provides  a 
correct  microscopic  description  of  the  nonlinear  response  of  a  molecule. 

Our  studies  of  quinoid  ring  structure  DCNQI  (sec  previous  section  of  Research 
Highlights)  motivated  us  to  perform  a  complete  theoretical  analysis  of  the  second  order 
properties  of  the  related  extended  ring  structure  DPHQ  [2-(4'-dicyanomethylene  bis 
cyclohexa-2,4,5^',5'-pcntacnylidine)-imidazolidine] 

)-o-o-o 

The  calculated  gas  phase  dispersion  curve  for  of  DPHQ  is  given  in  Figure  2.3.2.  It  is 
seen  that  the  magnitude  of  the  Px  value  is  unusuaUy  large,  achieving  nonresonant 

magnitudes  greater  than  1000x10'  esu  which  is  the  largest  value  ever  observed.  The 
origin  of  the  magnitude,  sign,  and  dispersion  of  resides  in  the  electron  density 

redistribution  accompanying  second  order  virtual  processes  in  the  dominant  excited  state  I. 
The  characteristic  features  of  this  principal  excited  state  I  are  illustrated  in  the  P|-Pg 
contour  diagram  of  Figure  2.3.3a,  along  with  the  corresponding  PfPg  contour  diagram 
for  MFC  shown  in  Figure  2.3.3b  for  comparison.  The  contour  diagram  illustrates  the 
highly  separated,  localized  nature  of  the  virtual  electron  density  distribution  of  DPHQ  along 
the  dipolar  x-direction  compared  to  the  modulated  distribution  in  MPC.  The  net  result  is 
that  although  the  molecular  length  along  the  dipolar  axis  is  roughly  the  same  for  MPC  and 
DPHQ,  the  corresponding  magnitudes  of  the  dipole  moment  Ap,*  are  quite  different,  equal 
to  4.4D  and  7.7D  respectively.  Similar  results  are  found  for  the  transition  moments  pj'g. 


12 


The  highly  correlated  electron  density  redistribution  is  the  primary  reason  for  the  unusually 
large  fijj  value  for  DPHQ. 

In  summary,  we  have  demonstrated  by  both  theoretical  and  experimental  studies  of 
extended  cyclic  ring  structures  that  our  theoretical  procedures  can  successfully  complement 
nonlinear  optical  experiments.  Moreover,  our  microscopic  description  for  in  terms  of 

highly  correlated  excitation  processes  is  essentially  correct  and,  thereby,  can  guide  I’he 
molecular  designs  of  new  structures  having  values  greater  than  1000x10  esu. 
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rv.  NONLINEAR  OPTICAL  PROCESSES  IN  CONJUGATED  POLYMERS: 

LARGE  NONRESONANT  THIRD  ORDER  ELECTRONIC  RESPONSES  IN  HIGH 
PERFORMANCE  UQUID  CRYSTAL  POLYMER  STRUCTURES 

Organic  and  polymer  structures  exhibit  unusually  large,  ultrafast  second  and  third 
order  nonlinear  optical  properties  in  a  large  number  of  structures,  phases,  and  states. 

Rapid  advances  in  the  field  can  be  achieved  through  further  development  of  stable,  high 
performance  polymer  structures  having  outstanding  secondary  properties  as  demonstrated, 
by  example,  by  high  performance  liquid  crystal  PBZ  polymers  (Figure  2.4.1).  For  two 
such  rigid  rod  polymer  systems,  PBI  and  PBT  (Figure  2.4.2),  our  third  harmonic 
generation  measurements  show  that  they  possess  large  non-iesonant  third  order  optical 
susceptibilities  whose  origin  resides  in  ultrafast,  lossless  excitations  of  highly  charge 
correlated  n-electron  states. 

Since  each  rigid  rod  polymer  chain  is  essentially  a  densely  packed,  linear  array  of 
linked  monomer  units,  a  repeat  unit  can  be  viewed  as  the  finite  chain  limit  to  the 
corresponding  infinite  polymer  chain.  Accordingly,  we  lepOTt  experimental  and  theoretical 
smdies  of  the  microscopic  third  order  susceptibilities  Yij)^|{-ci)4:©3, ©2,0)1)  of 

centrosymmetric  TBOZ  and  TBTZ  repeat  unit  structures  of  PBZ  polymers. 

The  nonresonant  macroscopic  third  harmonic  susceptibility  x^l^(-3©;©,©,©)  was 

determined  in  the  isotropic  liquid  phase  by  the  wedge  Maker  fringe  method.  For  an 
isotropic  liquid,  is  related  to  the  effective  microscpic  third  order  susceptibility 

Yl(-3©;©,©,©)  by 

2+nf,  „  2+nA 

x< -  N  ( -3-^3  (  PAl) 

where  N  is  the  number  density,  Yl  the  orientational  average  of  the  microscopic 
susceptibility  Yjjw(-3©:©,©,©)  and  the  refractive  index  terms  are  local  field  factors. 
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The  frequency  dependent  microscopic  7jj^j(-3©:ti),a),co)  shown  in  Figure  2.4.3 
was  theoretically  determined  by  SCF-MO-SDCI  calculations  including  singly  and  doubly 
excited  configurations  to  account  for  electron  correlations.  At  X,  ■  1 .907  p,m  (0.65  eV), 
the  calculated  Y|_  value  1 1  x  10'  esu  agrees  well  with  the  experimental  value  of  13  x  10' 
esu.  The  major  finding  is  that  the  microscopic  origin  of  YjjKi('3o);a),a),co)  is  determined 
by  virtual  excitations  to  highly  correlated  two  photon  rc-elcctron  states  through  the  first 
principal  one  photon  excited  state.  Important  features  of  the  highly  correlated  virtual 
processes  are  compactly  described  by  electron  density  correlation  diagrams  Pp^..  Figure 
2.4.4a  shows  Ppp.  for  the  ^A^  ground  and  excited  states,  and  Figure  2.4.4b  Pnn-  for 
the  ^B^  and  two  photon  ^A^  excited  states.  Thus,  we  also  demonstrate  the  necessity  of 
including  doubly  excited  configurations  in  third  order  calculations  as  well  as  the  important 
roles  of  symmetry  and  two  photon  states  in  third  order  processes. 
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Figure  2.4.3  Calculated  frequency  dependent  T|_(-3tt;oB,os.o).  Insert 

Experimental  and  calculated  optical  absorption  spectra.  Solid  Lme 
Theoretical:  dotted  line:  experimental  gas  phase;  dashed  line;  experimer‘a< 
solution  phase  in  dioxane. 
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Figure  2.4.4  Contour  diagram  of  transition  density  matrix 
ground  and  ^  By  excited  states  and  (b)  ^  By  one  photon  ind  ^  Ag 
excited  states  of  repeat  unit. 
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V.  NONLINEAR  OPTICS  AND  ORIENTATIONAL  DISTRIBUTION  FUNCTIONS 

OF  NEMATIC  UQUID  CRYSTAL  PHASES 


Relationships  between  the  optical  third  harmonic  susceptibility  ;5j^^|(-3to;cD,(0,o)) 
and  the  nonvanishing  order  parameters  <P2>  and  <P4>  for  nematic  liquid  crystals,  liquid 
crystal  polymers,  and  uniaxially  aligned  polymers  have  been  developed 


.  N  (8  +  25  <P2>  +  8n  <P4>) 

(2.5.1) 

.  N  (8  -  5  <P2>  +  3ii  <P4>) 

(2.5.2) 

iS  -  NS 

(2.5.3) 

where  8,  ^  and  q  are  explicit  summations  of  tensor  components  of  the  microscopic 
susceptibility  7ij|^,(*3(o;(o,co,to) . 

Third  harmonic  generation  measurements  (Figures  2.5.1  and  2.5.2)  were 
performed  on  the  model  system  nematic  N-(p-methoxybenzylidenc)-p-butylaniline 
(MBBA)  using  the  wedge  cell  technique  as  a  function  of  polarization  and  temperature  to 
demonstrate  the  simultaneous  determination  of  <P2>  and  <P4>  (Figures  2.5.3  and  2.5.4). 

The  results  are  in  agreement  with  previous  experimental  results  from  nuclear  magnetic 
resonance,  diamagnetic  anisotropy,  refractive  index,  and  Raman  scattering  studies.  Thus, 
we  have  successfully  developed  a  new  theoretical  formalism  and  experimental  methodology 
for  understanding  and  determining  the  relationships  between  third  order  nonlinear  optical 
processes  and  the  orientational  distribution  functions  for  nematic  liquid  crystals,  liquid 
crystal  polymers,  and  uniaxially  aligned  polymers.  Moreover,  the  formalism  and 
methodology  are  generalized  to  large  classes  of  partially  ordered  phases  and  to  second 
order  nonlinear  optical  processes. 
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VI.  SECOND  HARMONIC  GENERATION  AT  MILLIWATT  INCIDENT  POWERS 
IN  ORGANIC  STRUCTURES 

Theoietical  and  experimental  studies  have  demonstrated  the  important  features  of 
highly  charge  correlated  n-electron  states  in  organic  and  polymer  structures  that  comprise 
the  microscq)ic  origin  of  excq>tional  second  order  nonlinear  optical  responses 
5Cifk  •^)'  particularly  as  exhibited  in  second  harmonic  generation  (SHG)  and 

linear  electrooptic  (LEO)  properties.  As  an  important  example,  we  had  earlier  shown  that 
single  crystals  of  MNA  (2-mcthyl-4-nitroaniline)  possess  SHG  and  LEO  figures  of  merit 
SO  times  those  of  inorganic  dielectric  insulators  such  as  KDP  and  200  times  those  of 
semiconductors  such  as  GaAs  and  that  these  second  order  responses  are  primarily  purely 
electronic  at  dc  and  optical  frequencies. 

More  recently  we  provided  detailed  theoetical  and  experimental  analyses  of  the 
frequency  dependent  microscopic  second  order  response  3j||^(>2oD;co,Q>)  and  the  charge 

correlated  ic-electron  states  of  MNA.  Based  on  these  results,  we  now  report  theoretical 
and  experimental  results  of  visually  observable  phase  matched  (PM)  SHG  at  low  incident 
power  from  MNA  single  crystals.  PM  SHG  measurements  at  milliwatt  incident  power 
were  performed  on  x-y  and  x-z  cut  MNA  crystals  where  the  x-direction  optimizes  the 
microscopic  on-site  response  from  the  charge  correlated  p-electron  states.  For  example,  for 
a  60  mm  sample  interaction  length  in  the  x-y  configuration,  large  PM  SHG  is  easily 
observed  at  0.532  ^un  for  incident  Nd:YAG  (1.064  pm)  power  as  small  as  2  milliwatts. 

The  theoretically  calculated  PM  SHG  curves  for  Type  1  phase  matching  shown  in 
Figure  2.6.1a  were  confirmed  by  PM  SHG  experiments  at  normal  incidence  (Figure 
2.6.1b).  The  PM  SHG  curves  are  represented  on  a  grid  of  spherical  coordinates,  0  and  0, 
and  the  vertical  axis  is  the  effective  second  harmonic  component  There  is  a  maximum  at 
the  position  0  ■  50.1°,  0  «  1 .0°  that  was  extrapolated  to  the  position  0  »  51 .0°  and  o 
■  0°.  An  MNA  crystal  was  cut  and  pobshed  such  that  a  polarized  fundamental  laser  beam 
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at  X  ■  1 .064  (im  would  propagate  at  this  position.  The  observed  PM  SHG  locus  was 
followed  by  rotating  the  crystal  through  the  horizontal  and  vertical  axes. 

In  addition  to  fundamental  physical  understanding,  this  demonstration  of  readily 
observable  PM  SHG  at  low  incident  laser  power  also  serves  as  proof  of  principle  for 
applications  in  optical  signal  processing,  telecommunications,  and  data  storage  systems  that 
standardly  utilize  low  power  laser  diodes  as  optical  sources. 


Vn.  TOTAL  INTERNAL  REFLECTION  SPATIAL  LIGHT  MODULATOR  USING 
ORGANIC  ELECTROOPTIC  CRYSTAL  MNA. 

Nonlinear  optical  organic  and  polymeric  materials,  with  their  exceptionally  large 
nonlinear  optical  susceptibilities,  are  superior  in  many  aspects  to  commonly  used  inorganic 
materials.  This  clearly  has  provided  an  exciting  prospect  for  their  application  in  optical 
devices.  To  give  an  initial  assessment  of  the  implementation  of  these  organic  materials  in 
real  devices,  we  have  performed  a  thorough  study  of  an  optical  signal  processing  device, 
namely,  the  total  internal  reflection  (TIR)  spatial  light  modulator  (SLM)  and  conxpare  the 
results  with  the  performance  characteristics  of  this  device  based  on  a  standard  inorganic 
material  as  the  active  medium.  The  comparison  was  made  between  crystal  MNA  (2- 
methyl-4-nitroaniline),  a  well  characterized  organic  electrooptic  and  nonlinear  optical 
material,  and  LiNb03,  a  widely  employed  inorganic  electrooptic  crystal. 

The  basic  structure  of  TIR  SLM  is  shown  in  Figure  2.7.1.  The  device  consists  of  a 
VLIC  integrated  circuit  having  serial-to-parallel  addressing  electronics  that  provide  data  to 
several  thousand  drive  transistors  which  then  connect  to  several  thousand  parallel 
electrodes  on  the  surface  of  the  chip.  This  chip  is  pressed  against  an  electrooptic  crystal  so 
that  fringing  fields  existing  between  electrodes  are  proximity  coupled  into  the  crystal 
(Figure  2.7.2).  These  fringing  fields  change  the  index  of  refraction  of  the  crystal  through 
the  linear  electtooptic  effect  To  interact  with  these  fringing  fields,  the  device  is  illuminated 
with  a  sheet  of  light  which  is  collimated  in  the  horizontal  axis  and  brought  to  a  sheet  focus 
on  the  surface.  As  the  light  approaches  this  surface  and  recedes  from  this  surface,  it  picks 
up  a  phase  shift  dependent  on  the  index  of  refraction  change  produced  by  the  fringing 
electric  fields.  By  using  a  standard  Schlieren  imaging  system,  this  phase  modulated  profile 
can  be  converted  into  an  optical  image  intensity  profile. 

Performance  of  the  TIR  SLM  is  evaluated  for  MNA  and  LiNbOj  using  reported 

electrooptic  coefficients.  Our  results  demonstrate  that  to  induce  the  same  values  of  optical 
phase  modulation,  the  required  driving  voltage  for  a  TIR  device  using  MNA  as  the 


3 


modulator  crystal  is  reduced  by  approximately  a  factor  of  2.5,  ccnnpared  with  that  using 
LiNbOj  as  modulator  crystal.  Correspondingly,  under  the  same  driving  voltage  pattern, 

the  diffraction  efficiency  will  be  enhanced  by  a  factor  of  more  than  six  for  a  TIR  device 
using  MNA  as  compared  with  that  using  LiNb03  as  naodulator  crystal.  These  results  for 

MNA  also  show  that  poled  electrooptic  polymers,  which  are  formed  basically  to  emulate 
MNA  crystals  but  in  polymer  form,  can  additionally  be  implemented  in  TIR  SLM 
configurations  with  figures  of  merit  comparable  to  and  greater  than  those  of  an  MNA  based 
SLM. 


32 


VIIL  SYMMETRY-CONTROLLED  ELECTRON  CORRELATION  MECHANISMS  OF 
THIRD  ORDER  NONLINEAR  OPTICAL  PROPERTIES  OF  CONJUGATED  LINEAR 
CHAINS 

Third  order  nonlinear  optical  polarizadons  created  by  interactions  between  intense 
laser  fields  in  optical  media  result  in  a  large  variety  of  self-action  phenomena  that  are 
important  to  fundamental  physics  and  materials  research  as  well  as  optical  device 
technologies.  Self-action  effects  are  the  basis  for  controlling  light  by  light  in  its 
propagation,  polarization,  amplitude,  and  frequency.  Important  examples  include  harmonic 
generation,  self-phase  modulation,  phase  conjugate  wave  generation,  and  processes  such 
as  optical  bistability  that  are  important  to  optical  signal  processing  and  data  storage. 

It  is  well  established  that  conjugated  organic  and  polymer  structures  exhibit 
unusually  large,  ultrafast  third  order  nonlinear  optical  properties.  These  properties  arise 
from  virtual  p-election  excitations  occurring  within  condensed  assemblies  of  individual 
molecular,  or  polymer  chain,  sites  which  act  as  macroscopic  sources  of  nonlinear  optical 
response.  However,  until  the  present  study,  no  suitable  fundamental  understanding  existed 
for  these  unusual  third  order  properties. 

This  year  through  a  series  of  experimental  and  theoretical  studies,  we  have 
successfully  developed  a  tmcroscopic  description  of  the  mechanism  and  a  corresponding 
set  of  physical  principles  for  understanding  third  order  virtual  excitation  processes  that 
occur  in  the  p-electron  states  of  conjugated  linear  chains  such  as  polyenes  and  polyenyres 
(Figure  2.8.1).  We  have  shown  that  symmetry  controlled  correlation  effects  determine  the 
properties  and  behavior  of  the  third  order  nonlinear  optical  susceptibilities  of  these 
structures.  The  highly  correlated  nature  of  the  third  order  virtual  excitations  for  the  third 
harmonic  susceptibility  diroctly  represented  by 

contour  diagrams  of  transition  density  matrices  p^p..  Contour  diagrams  for  Pp^,.  of  the 
ground  state  and  highly  electron  conelated  6' A^  state  with  the  dominant  first  one  photon 

1  state  are  shown  in  Figure  2.8.1  where  solid  and  dashed  lines  correspond  to  opposite 
signs  of  Ppp..  The  contour  diagram  of  Pgj  'b  the  g  1  transition  results  in  a 
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modulated  charge  redistribution  and,  correspondingly,  a  relatively  small  contribution  to 
Yxxxx-  ^  sharp  contrast,  ^10  for  the  virtual  transition  between  the  and  l^B^^ 

states  produces  a  large  charge  separation  along  the  x-chain  axis  direction  and  an  associated 
large  transition  moment  of  13.2D  that  dominates  the  contributing  term  to  Yxxxx- 

The  calculated  dispersion  curve  for  y^^^j((-3tD;o),co,to)  for  a  principle  example  of 

octatetraene  (OT)  having  eight  carbon  atom  sites  (N)  is  shown  in  Figure  2.8.1.  Over  the 
near  infrared  and  visible  ranges,  a  first  peak  at  1.47  eV  (X  ■  0.84  pm)  appears  due  to 
the  3(0  resonance  of  the  l^B^^  state,  and  a  second  peak  at  2.08  eV  (X  ■  0.60  pm)  is  from 
the  2^Ag  state  which  is  significant  on-resonance  even  though  it  makes  a  very  small 
contribution  off-resonance. 

The  general  features  of  the  origin  of  Yj|^|{-3(D:a),(o,co)  for  OT  just  described  are 

common  to  all  the  chain  lengths  from  N  =  4  to  16.  Importantly,  there  is  a  unifying  general 
result  that  at  fixed  frequency,  Yxxxx  differing  chain  lengths  exhibits  a  dramatic  but 
well-defined  increase  with  increased  chain  length.  The  theoretical  values  for  it  a 
fundamental  photon  energy  of  0.65  eV  (X  »  1 .907  pm)  are  plotted  against  the  number 
of  carbon  atom  sites  N  on  a  log-log  scale  in  Figure  2.8.2  and  exhibit  a  power  law 
dependence  upon  N  with  an  exponent  of  5.4  which  results  from  three  length-dependent 
factors.  With  increased  chain  length,  there  are  lower  excitation  energies,  larger  transition 
dipole  moments,  and  increased  number  of  contributing  excited  state  terms. 

In  summary,  we  have  developed  a  microscopic  mechaiusm  for  third  order  nonlinear 
optical  properties  of  conjugated  linear  chains  that  allow  detailed  identification  of  the 
significant  contributions  to  any  third  order  process.  Such  a  microscopic  description  further 
allows  development  of  guidelines  for  the  molecular  design  of  new  organic  nonlinear  optical 
materials. 


MICROSCOPIC  ORIGIN  OF  r(-3u;u,u.ut 
FOR  LINEAR  CONJUGATED  CHAINS 
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IX.  fflGH  RESOLUTION  STRUCTURAL  STUDIES  OF  ULTRATHIN 

MULTILAYER  FILMS  CONTAINING  A  DISUBSTITUTED  DIACETYLENE 
POLYMER 


Disubstituted  diacetylenes  can  be  polymerized  in  the  solid  state  to  form  long-chain 
polymers.  The  extended  7i-electron  system  of  the  conjugated  polymer  is  highly 
polarizable,  and  the  excited  tt-electron  states  give  rise  to  extremely  large  microscopic 
nonlinear  q}tical  susceptibilities.  These  microscopic  susceptibilities  can  result  in  large 
macroscopic  nonlinear  optical  properties,  which  mediate  the  observed  nonlinear  optical 
processes  of  the  material  medium.  Recently,  interest  has  centered  on  nonlinear  optical 
properties  of  Langmuir-Blodgett  multilayer  films  of  amphiphilic  disubstituted  diacetylenes 
because  such  films  provide  ordered  polymer  structures  in  ultrathin  film  forms. 

We  have  recently  completed  high-resolution  x-ray  scattering  structural  smdies  on 
both  periodic  and  nonperiodic  ultrathin  multilayer  films  containing  from  one  to  five 
bilayers,  each  bilayer  of  either  arachidic  acid  or  pentacosa-10,12-diynoic  acid  (a 
disubstituted  diacetylene  monomer).  The  monomer  was  UV-polymerized  in  the  films  into 
either  its  blue  form  "DAPb"  or  its  red  form  "DAPr".  Some  of  the  results  are  summarized  in 
Figure  2.9.1  below. 

The  figure  presents  a  schematic  view  of  our  current  knowledge  of  the  structure  of 
periodic,  ultrathin  multilayer  films  containing  a  disubstituted  (poly)  diacetylene.  The 
double  lines  running  parallel  to  the  plane  of  the  multilayer  are  intended  to  represent  the 
backbone  of  the  extended  polyene  chain,  which  links  adjacent  monomeric  units  together. 
The  first  deposited  monolayer  is  depicted  as  being  composed  of  well-ordered  all-trans 
hydrocarbon  chains  with  their  methyl  ends  distributed  over  a  small  distance  5^  along  the  z- 

axis.  Once  a  bilayer  has  been  deposited,  the  underlying  OTS  layer  is  shown  as  disordered 
via  chain  configurations  containing  "kinks"  and  "jogs"  in  such  a  way  that  the  variable  Cjg 

chain  lengths  accommodate  the  "rough"  class  surface.  The  average  tilt  angle  t  of  the 
hydrocarbon  chain  axis,  relative  to  the  normal  to  the  layer  plane  is  _  35°  for  DAPb  and 

for  DAPr.  Each  interiOT  bilayer  for  the  DAPb  multilayer  is  shown  as  containing  well- 
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ordered,  essentially  all-trans  chains  with  a  small  spread  of  the  position  of  the  methyl  end 
groups  along  the  z-axis  of  5i .  Asa  result,  the  interlayer  ordering  in  both  periodic  and 

nonperiodic  multilayer  films  containing  DAPb  is  long-range.  The  interior  bilayers  of  the 
DAPr  multilayer  are  shown  as  having  well-ordered  all-trans  carboxyl  end  side-chains, 
while  the  methyl  end  side-chains  are  disordered  via  "kinks''  and  "jogs"  with  their  methyl 
end  groups  spread  over  a  larger  length  52-  As  a  result,  the  interlayer  ordering  in  both 

periodic  and  nonperiodic  multilayer  films  containing  DAPr  is  only  short-range.  The  last 
monolayer  is  shown  with  significant  disorder  of  the  hydrocarbon  chains  distributing  the 
methyl  end  groups  over  a  much  larger  length  e  along  the  z-axis.  The  double  lines 
representing  the  polyene  backbone  are  shown  as  dashed  lines  for  the  last  monolayer  as  the 
degree  of  disorder  may  reduce  the  extent  of  polymerization. 
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Figure  2.9.1 


X .  FREQUENCY  DEPENDENCE  OF  THE  OPTICAL  KERR  EFFECT  AND  THIRD 
ORDER  ELECTRONIC  NONLINEAR  OPTICAL  PROCESSES  OF  ORGANIC 
SYSTEMS. 

Self-action  effects  in  third  order  nonlinear  optical  processes  provide  the  basis  for 
controlling  light  with  light  and  are,  thereby,  important  to  both  fundamental  physics  and 
optical  device  technologies.  Third  order  nonlinear  optical  processes  include  third  harmonic 
generation,  optical  Kerr  effect,  phase  conjugate  wave  generation,  and  optical  bistability. 

In  material  media,  mechanisms  involving  electronic,  lattice,  and  thermal  excitations 
all  contribute  in  experimental  determinations  of  macroscopic  These  contributions 

occur  over  a  broad  frequency  range  extending  from  the  near  infrared  to  the  ultraviolet 
Highly  precise  experimental  methods  are  therefore  required  for  accurately  determining  each 
term  possibly  contributing  in  a  experiment  Third  harmonic  generation  (THG)  for 

example,  measures  only  the  electronic  term  and  is  severely  limited  in  frequency  range 
because  either  the  fundamental,  the  generated  wave,  or  both  often  are  absorbed  by  nearby 
one  or  two  photon  resonances.  Degenerate  four  wave  mixing  (DFWM)  configurations  are 
fixed,  single  frequency  measurements,  and  the  output  signal  includes  all  contributing  terms 
that  are  not  resolved  unambiguously  even  down  to  subpicosecond  tunescales.  Thus, 

values  determined  by  THG  or  DFWM  methods  are  severely  limited  and  oftentimes 
unreliable. 

We  have  developed  new  experimental  methods  based  on  the  frequency  dependence 
of  the  optical  Kerr  effect  and  nonlinear  optical  interferometry  for  accurately  determining  the 
magnitude,\sign,  and  dispersion  of  This  year  we  have  completed  an  extensive 

series  of  experimental  and  theoretical  studies  of  the  frequency  dependent  optical  Kerr  effect 
that  enables  frequency  dependent  j^i^Vesults  to  be  accurately  obtained  from  optimized 

pump-probe  techniques. 

The  optical  Kerr  effect  (OKE)  is  a  four  wave  mixing  process  in  which  an  intense 
optical  pump  beam  of  frequency  co.|  induces  depolarization  in  a  probe  beam  of  frequency 
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0)2  because  of  the  thiixl  order  nonlinear  optical  polarization,  (0)2).  In  an  isotropic 
medium,  (0)2)  is  expressed  as 


(©2)  »  Xjjk/(-0)2:0)i.-0)i.0)2) 


,(3.) 


(2.10.1) 


where  the  third  order  nonlinear  optical  susceptibility  and  and  are  the 

pump  and  probe  fields  respectively.  At  frequencies  well  above  molecular  vibration  and 
rotation  modes  but  below  electronic  resonances,  (-©gio)^  ,-q)^  ,0)2)  has 

contributions  due  to  both  electronic  excitations  and  lattice  displacements  which  in  the  Bom- 
Oppenheimer  approximation  is  expressed  as  the  siun  of  the  quantum  averages  Gjji^i  and 


XiW  -  <^ijkl  +  dijkl  (2.10.2) 

where  and  dy^  are  the  third  order  electronic  and  lattice  contributions,  respectively; 
both  terms  are  assumed  to  be  fast  relative  to  experimental  time  scales.  In  such  a  frequency 
regime,  because  of  nearby  electronic  resonances  at  higher  frequency,  the  dominant 
fiequency-dependent  term  is  the  ultrafast  electronic  contribution,  and,  hence,  fiequency- 
dependent  measurements  of  the  total  OKE  response  allow  direct  determination  of  the 
electronic  and  lattice  contributions  to  fro®  OKE  data  alone. 

In  our  experimental  studies,  the  OKE  coefficient  is  measured  utilizing  a  pump- 
probe  configuration  with  the  pump  and  probe  beam  polarization  oriented  at  45°  with  respect 
to  each  other  (see  Figure  2.10.1).  The  pump  beam  induces  a  refractive  index  anisotropy  in 
the  Kerr  medium,  thus  resulting  in  a  rotated  elliptical  polarization  for  the  probe  beam.  The 
polarization  change  is  detected  using  an  appropriate  arrangement  of  cross  polarizers.  In  the 
experiment,  the  pump  beam  employed  is  a  pulsed  Nd:YAG  laser  (1.064  pm)  while  the 


piobe  beam  is  an  argon  i(Mi  laser  or  a  cw  ring  dye  laser  to  permit  tunability  and  thus 
facilitating  a  dispersion  measuiemenL  An  optical  gate  system  was  designed  to  reduce  the 
background  noise  to  the  PMT  detector  system  due  to  the  cw  source. 

The  frequency  dispersion  of  the  OKE  susceptibility  a  -Q)2;co, ,-co^,a)2)  was 

measured  for  a  number  of  systems.  Results  for  the  model  Kerr  liquids  benzene  and 
nitrobenzene  are  shown  in  Rgure  2.10.2  .  In  general,  the  >^2) 

composed  of  a  resonant  term  and  a  nonresonant  term.  The  resonant  term  is  due  to  the  two- 
photon  resonance  of  the  sum  (o).|  +CO2)  or  difference  -tD2)  frequency  of  the  pump 
(co^)  and  probe  {<02)  beam,  while  the  nonresonant  term  is  in  part  due  to  the  molecular 

reorientation  contribudon.  We  have  demonstrated  for  the  frrst  time  that  since  the 
experiment  measures  the  absolute  value  of  the  ;{|ki^‘0>2;ci>^  ,~o)^  ,0)2),  interference 

between  the  resonant  and  the  noniesonant  terms  gives  rise  to  a  characteristic  Fano-like 
resonance  peak,  as  can  be  seen  in  the  experimental  data  of  Figure  2. 10.2.  From 
comparison  with  gas  phase  linear  absorption  data  and  theoretical  calculations  of 

>^2)  observed  dispersion  data  can  be  thoroughly  understood  and 
the  resonant  electronic  contribution  determined.  OKE  measurements  have  the  principal 
advantage  that  since  both  the  pump  and  the  probe  frequencies  are  far  off  resonance, 
negligible  single  photon  absorption  occurs  and  there  is  little  or  no  depletion  of  either  beam. 
This  allows  a  thorough  study  of  the  resoiuuit  electronic  contributions  to 


In  summary,  an  experimental  and  theoretical  study  of  OKE  dispersion  in  isotropic 
media  has  been  developed.  A  theoretical  analysis  based  on  Maxwell's  equations  provides 
the  connection  between  experimentally  measured  quantities  and  1 1  .^^2)  • 

An  optimized  pump-probe  depolarization  technqiue  was  developed,  and  conventional 
theoretical  approximations  thoroughly  studied  The  OKE  dispersion  data  is  used  to 

determine  the  electronic  and  lattice  contributions  to  Xj*ki^-W2’“i’*®i’®2) 

frequency  electronic  resonances  and  can  show  the  relative  strength  of  the  electronic  term, 


Ojji^i,  from  different  excited  states.  The  dispersion  results  for  7t-electron  systems  with 
electronic  transitions  probed  in  the  presence  of  a  continuum  background  provided  the  first 
observation  of  a  characteristic  Fano-like  peak  due  to  interference  between  resonant  and 
nonresonant  terms.  The  general  formalism  presented  in  this  study  can  be  adopted  for  the 
study  of  resonant  and  nonresonant  processes  in  any  other  media  of  interest;  and,  the 

experimental  data  can  be  compared  to  microscopic  description  based  upon  many-electron 
theory  and  calculations. 
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Experimental  layout  where  HWPi  ia  a  half  wave  plate  at  i.06Mm,  Pi 
and  P2  polarizera.  Ai  and  a2  anaiizera.  BSi  and  BS2  beamapiittera,  MOMi  and 
MOM2  a  motor  driven  mirror,  KOP  a  KOP  mixing  cryatal.  PO  a  photodiode,  F  i 
and  P2  flltera  and  PMT’a  are  photomultiptier  tubea. 


neoaacoNO  oxa  ixnaiwiiT 


Figure  2.10.1 
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Oiaptrtion  data  for  nitrobanzaaa  ia  givan  by  tha  solid  points  (scaia 
on  latt  ordinata)  and  banzana  data  ia  givan  by  tha  opan  points  (scaia  on  right 
ordinata). 


Figure  2.10.2 


XI.  TUNABLE  ODHERENT  LASER  SOURCE  FROM  STIMULATED 
ROTATIONAL  RAMAN  SCATTERING 

We  have  coaspleted  a  detailed  study  of  stimulated  rotational  Raman  scattering 
(SRRS)  using  a  Nd:YAG  laser  and  demonstrated  that  it  can  be  used  as  a  convenient,  low 
cost  coherent  light  source  from  0.4  |im  to  2.2  |jm.  Stimulated  Raman  scattering  (SRS)  is 
widely  used  to  frequency  shift  the  1.064  |im  output  from  the  Nd:YAG  laser  using,  for 
example,  hydrogen  gas  where  the  vibrational  excited  state  provides  an  energy  shift  of  4155 
cm'^  that  results  in  four  additional  coherent  sources  in  the  visible  and  near  infi^ued  regions 
(1.907, 0.737, 0.564  and  0.457  |im).  The  rotational  Raman  shift,  being  much  smaller  in 
energy,  results  in  more  closely  spaced  lines.  The  rotational  level  spacing  AE  of  a  diatomic 
molecule  can  be  obtained  by  considering  a  quantum  vibrating  rotator,  and  is  given 
approximately  by  the  expression 

AE  -  [Be-ao(v+1/2)]J(J+1)  -  [De-Pe(''+1/2)]j2(J+1)2  (2.11.1) 

where  J  is  the  total  angular  momentum,  v  the  oscillator  quantum  number,  and  Bg,  oc^,  Dg 
and  3g  are  constants  for  the  molecule.  With  the  selection  mle  AJ=2,  the  energy  spacing 
for  the  rotaticmai  levels  is  586.9  cm'^.  This  effect,  in  combination  with  the  vibrational 
Raman  lines,  produces  closely  spaced  discrete  lines  spanning  from  0.4  mm  in  the  visible  to 
2.2  ^m  in  the  near  infrared. 

A  high  power  pulsed  NdcYAG  laser  was  used  to  pump  a  high  pressure  H2  gas  cell. 
With  the  appropriate  input  polarization  up  to  five  Stoke  lines  can  be  easily  observed  using  a 
moderate  input  power.  This  effect,  together  with  the  vibrational  Raman  lines  produces 
more  than  20  discrete,  well-defined  new  frequency  lines  in  the  visible-near  infiared  region 
using  only  a  single  frequency  pump  source.  In  one  demonstration,  we  have  utilized  these 
SRRS  lines  in  measurements  of  the  dispersion  of  second  and  third  harmonic  generation  in 
organic  systems.  Results  of  experimental  studies  of  the  power  dependence  of  the  SRRS 
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lines  at  various  frequencies  on  gas  pressure,  polarization,  and  power  of  the  input  source  are 
shown  in  Figures  2.1 1.1  to  2.1 1.3.  The  simplicity  and  reliability  of  SRRS  makes  it  a 
convenient,  low  cost,  tunable  source  of  visible  and  near  infrared  coherent  radiation.  The 
case  of  operation,  for  instance,  may  be  important  to  issues  of  sensor  protection. 
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Figure  2.11.1 
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Figure  2.11.2 
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POWER  DEPENDENCE 


Figure  2.11.3 
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Ill 


EXTENDED  DESCRIPTION 


I .  NONLINEAR  OPTICAL  PROCESSES  IN  LOWER  DIMENSIONAL 
CONJUGATED  STRUCTURES 

1.  INTRODUCTION 

Ccmjugated  n-electron  organic  and  polymer  structures  are  now  well-known  to  exhibit 
unusually  large  noniesonant  macroscopic  second  order  and  third  order 

X^^(-(n4;ci)i,Ci)2>ci>3)  nonlinear  t^tical  susceptibilities,^^^  and  their  microscopic  origin  and 

mechanism  can  be  successfully  described  by  quantum  field  theory  of  many-electron  systems  in 
one  and  two  dimensions.  In  this  description,  as  the  spatial  dimensionality  of  the  many- 

electron  system  is  effectively  lowered,  the  motion  among  the  many-electrons  becomes  highly 
coiielated.^^^^  Electron  correlation  plays  a  major  role  in  determining  Xj^^  and  Xj^  in 

ctmjugated  ic-electron  systems,  and,  as  we  shall  review,  its  effects  cannot  be  neglected  in 
properly  accounting  for  the  nonlinear  optical  properties  of  conjugated  structures. 

Reduced  dimensionaiity  and  quantum  confinement  occur  in  material  structures  having 
effective  confineoaent  lengtii  scales  L  of  less  than  one  to  several  tens  of  nanometers.  The 
opposite  limit  occurs  in  regular  three  dimensional  bulk  structures  where  L  >  100  nm,  and  the 
electrons  behave  as  a  Fermi  gas  of  weakly  interacting  particles.  In  this  case,  the  electron 
motions  are  usually  weakly  correlated  and  well  described  by  single-particle  theory  in  the 
effective  mass  approximation. 

Conjugated  x-electron  organic  and  polymer  structures  are  viewed  as  lower  dimensional 
systems  in  which  quantum  confinement  effects  are  important  For  instance,  in  conjugated 
linear  chains,  the  many-electron  system  is  confuted  in  two  dimensions  with  L  <  0.5  nm 
transverse  to  the  linear  chain  axis,  and  thus  the  x-electrons  are  delocalized  in  their  motion  oni> 
in  one  dimension  along  the  longitudinal  chain  axis.  Asa  result,  the  repulsive  Coulomb 
interactions  between  electrons  are  strong,  and  electron  motion  becomes  highly  correlated 
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wherein  the  motion  of  any  single  electron  depends  on  all  the  other  remaining  many  electrons. 
The  origin  of  the  nonresonant  Xyu  f<^  these  and  related  polymer  structures  then  appears  to 

reside  in  strong  correlation  behavior  in  virtual  two-photon  7C-electron  states. This 
microscopic  many-electron  coneladon  description  of  nonlinear  optical  properties  of  conjugated 
structures  is  described  in  section  2  of  this  paper.  We  first  review  the  nonlinear  optical 
properties  of  linear  polyenes  both  as  a  function  of  conformation  and  lengdi  as  prototype  one 
dimensional  cmjugated  chains.  The  nonresonant  microscopic  nonlinear  optical  susceptibility  is 
a  strong  function  of  the  length  of  the  polyene  chain  and  not  significantly  dependent  on 
structural  confonnadmi. 

The  role  of  dimensionality  in  nonlinear  optical  properties  is  further  examined  by  a 
theoretical  study  of  a  two  dimensional  conjugated  cyclic  structure  with  L  <  0.5  nm 
perpendicular  to  die  molecular  plane.  Through  comparison  with  the  linear  chain  results,  the 
microscopic  origin  of  the  nonlinear  optical  properties  is  shown  to  be  similar  for  the  two  cases. 
A  remarkable  effect  of  increased  dimensionality,  however,  is  the  decrease  of  the  isotropically 
averaged  third  order  susceptibility  due  to  an  actual  reduction  in  the  effective  length  available 
for  the  7C-electrons  to  respond  to  an  applied  optical  electric  field. 

In  section  3,  we  present  experimental  measurements  of  the  resonant  nonlinear  optical 
properties  of  thin  films  formed  of  large,  two  dimensional  conjugated  cyclic  structures. 
Saturable  absorption  is  observed  in  spin  coated  thin  films  of  a  naphthalocyanine  derivative 
having  an  exceptionally  large  Q-band  x-electron  absorption  (o  of  1  x  lO^  cm-')  centered  at 
810  nm.  No  unsaturable  absorption  background  is  observed.  For  laser  frequencies  near  810 
nm,  we  find  an  intensity  dependent  refractive  index  (02)  of  1  x  10"^  cm^/kW.  Based  on  the 


results  of  the  saturable  absoipdon  studies,  we  then  demonstrate  electronic  absorptive  opdcal 
bistability  in  a  thin  film  Fabty-Perot  etalon  at  nanosecond  time  scales. 


2 .  MICROSCOPIC  DESCRIPTION  OF  NONRESONANT  7ijkl(-oi4;o)i  ,(02,0)3)  FOR 

ONE  AND  TWO  DIMENSIONAL  STRUCTURES 


The  macroscopic  nonlinear  t^tical  properties  of  organic  molecular  and  polymer 
structures  in  condensed  states  are  best  described  starting  from  the  individual  responses  of 
isolated  molecular  or  polymer  chain  units.^^^  This  approach,  which  clarifies  the  origin  of 
and  is,  however,  not  i^)plicabie  to  other  material  classes,  such  as  inorganic 
semioonductors,  where  it  is  difficult  to  identify  analogous  microscopic  sources  of  nonlinear 

polarization.  For  the  moment  neglecting  intermolecular  intetactions,  if  tire  nonlinear 
susceptibility  of  an  isolated  molecule  is  known,  then  or  Xi^  of  the  macroscopic  ensemble 

of  molecules  is  determined  by  the  orientational  distribution  function  of  the  independent  units. 
Local  field  factors  must  also  be  included  to  account  for  the  effect  of  the  dielectric  environment 

on  the  electric  field  strengtii  at  the  molecular  site.  The  macroscopic  frequency-dependent 
second  and  third  order  susceptibilities  xj^(*<i>3;(i)i,(i>2)  and  Xiju('0>4;o)i, (02,(03)  can  then  be 

expressed  in  terms  ci  the  molecular  susceptibilities  Pijk(-(03;(0i,(D2)  and  Ytjki('(04:(0i, (02,(03) 
ast'31 


xS^(-<«3;o>,.a,2)  -  N„  Rj^.  Rj*.  R^V  f 


(1) 


and 
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(2) 


Xig(-<04;«)i,0):,(03)  -  N„  I  Ri^.  RjJ.  RkV  R,'. 

S*“  A 


where  Nu  is  the  number  of  unit  cells  per  unit  volume,  the  summation  is  over  all  molecules  in 
the  unit  cell,  R  is  a  rotation  matrix  describing  the  orientation  of  each  molecule  in  the  unit  cell, 
and  f  is  the  frequency  dependent  local  field  factor.  The  description  of  the  macroscopic 
nonlinear  optical  response  is  thus  reduced  to  an  understanding  of  the  microscopic  second  and 
third  order  susceptibilities  3ijk  and  Tiju  and  knowledge  of  the  orientational  distribution  of  the 
molecular  units  in  the  condensed  phase.  As  special  cases,  isotn^ic  gases  and  liquids  reduce 
Eqs.  (1)  and  (2)  to  simpler  forms.  For  example,  for  an  isotropic  ensemble 

Xnn(*<04;0)i,0)2,0>3)  *  N  f“^  f“^  f^  Yg(-<04;ci>i,G)2,tU3)  (3) 

where  N  is  the  number  density  of  molecules  and  yg  is  the  isotropically  averaged  susceptibility 
defined  by 


Yg  *  t[  J  Yiiii  T  ’ 

•^1  i#3 

where  the  indices  i  and  j  represent  the  Cartesian  coordinates  x,  y,  and  z. 

The  above  formalism  is  strictly  appropriate  only  for  gases  such  that  the  mean 
intermolecular  distance  is  large  enough  for  the  interaction  between  molecules  to  be  negligible 
In  condensed  states,  some  modifications  are  required  to  include  the  effects  of  intermolecular 
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interactions.  In  order  to  account  for  residual  interactions  with  neighboring  molecules, 
|5ijk(*0)3;ci)i,CD2)  and7ijki(-0)4;(0i,(i>2><i>3)  can  be  "dressed"  to  extend  their  applicability  beyond 
isolated  molecular  units.  For  instance,  in  dc-induced  second  harmonic  generation  (DCSHG) 
studies  of  dilute  soludons,^^^*^^^  it  was  found  that  inclusion  of  dipole-mediated  interactions 
between  solvent  and  solute  molecules  provided  excellent  agreement  of  die  theoretical  values 
for  an  otherwise  independent  molecule  with  experimental  measurements  made  in  solution.  In 
a  more  recent  smdy^^^^’^^^  of  phase-matched  second  harmonic  generation  in  single  crystals  of 
2-methyl-4-nitroaniline  (MNA),  the  existence  of  closely  oriented  pairs  of  molecules  in  the  unit 
cell  crystal  structure  required  treatment  of  an  MNA-MNA  pair  as  the  fundamental  source  of 
second  order  response.  Because  of  the  strong  mutual  effect  of  each  molecule  on  the  electronic 
structure  of  the  odier,  calculation  of  ^jk(-2a);(D,(o)  for  an  appropriately  oriented  pair  of  MNA 
molecules  was  necessary  in  order  to  obtain  agreement  with  experiment  In  any  case,  the  fact 
that  the  intranaolecular  interaction  energy  is  much  stronger  than  the  intermolcular  interaction 
energy  in  organic  molecular  crystals,  liquids,  solutions,  and  polymer  thin  films  means  that  the 
starting  point  for  understanding  die  macroscopic  nonlinear  optical  properties  of  diese  various 
condensed  phases  lies  in  an  accurate  description  of  the  microscopic  response  fiom  an  isolated 
molecular  unit 

The  general  thecnetical  expression  for  the  components  of  the  microscopic  third  order 
susceptibility  tensor  7iju(*ca4;a)i, 012,(1)3)  is  derived  from  time-dependent  quannim 
electrodynamic  perturbation  theory.  In  order  to  avoid  secular  divergences  that  would  occur 
when  any  subset  of  the  input  frequencies  sums  to  zero,  one  employs  the  Bogoliubov- 
Mitropolsky  method  of  averages.^^®^’^*^^  Owing  to  dispersive  effects,  yijki(*(04;0)i,ci)2,(‘>3)  is 
dependent  on  the  input  and  output  frequencies  involved  for  each  of  the  various  possible 
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nonlinear  optical  phenomena.  For  the  particular  case  of  third  harmonic  generation,  for 
example,  one  obtains 


4K3  nin2n3  (0)n3g-3(0)(a)n2g-2t0)(0)n,g-Ci)) 


^  Pjkl[r|n3  tfijiia  r^^m  filtgl 

(0)n3g+0))(0)n2g-2(D)(a)n,g-a>) 


Pjkl[r^n3  t‘n2ni  ^nig] 

(Q)n3g+<i>)(t0n2g+20))(Ci)n,g-C0) 


^  Pjkl[rin,  r^3n2  fn^ni  j 

(o>n3g+o))(con2g+2(0)(a)nig+3<o) 


I' 

nin2 


j  Pjkiifgnz  ^iiis  fgiii  t‘n\g] 

(0i)n2g-3tD)(0>n2r®)(“nir“) 


^  Pjklfrjna  Tnag  tgni  rn\g]  ^  Pjkltrjna  rfeg  r^,  ftiig] 
(0)n2g’0>)(<Dn,g+<D)((0n,g'0))  (tDn2g+3(i»(Oln2g-KD)(a)B,g-KD) 


^  Pjkl[f^n2  fflag  ^gni  ^n\g3  j  j 
(a)h2g-Ko)(a)n,g-o)Xti)|„g-««» 


(5) 


where  rn|n2  **  tnatrix  element  <nil  r*  In2>  ( r'  =  r*  -  <r‘>gg),  ITtOng  is  the  excitation  energy 
of  state  n,  the  prime  on  die  summations  indicates  the  ground  state  is  omitted,  and  Pju  denotes 
the  sum  over  all  permutations  of  those  three  indices.  We  will  be  concerned  here  with  the  case 
where  all  of  the  optical  frequencies  are  above  the  molecular  vibrational  and  rotational  modes 
but  below  the  electronic  excitation  energies  so  that  the  nonlinear  optical  response  is  strictly 
electronic  in  origin.  In  addition,  for  conjugated  organic  structures,  y  is  dominated  by  the 
delocalized  tc-electron  contributions  which  in  general  have  both  larger  transition  dipole 
tiximents  and  lower  transition  energies  than  the  o-clectron  excitations.  Thus,  with  accurate 


description  of  both  the  excitation  energies  and  transition  moments  of  the  x-electron  manifold, 
one  can  calculate  the  firequency  dependence  of  each  of  the  different  third  order  nonlinear 
optical  processes  using  expressions  such  as  Eq.  (S). 

The  theOTedcal  method  employed  to  achieve  proper  description  of  the  7C-electron 
manifold  of  conjugated  organic  naolecular  structures  for  calculation  of  7ijki('0)4;(i)l,0)2,C03)  has 
been  described  previously.^^^^  It  consists  of  a  multiply-excited  configuration  interaction 
calculation  (SCF-MO-SDCI)  performed  on  a  molecular  orbital  basis  obtained  through  self- 
consistent  field  theoy  in  the  rigid  lattice  CNDO/S  approximation.  This  method,  which  has 
previously  been  proven  successful  in  calculations  of  Pijk(-Q>3;ti>i,n)2)  for  a  variety  of 
conjugated  structures  will  not  be  discussed  here  other  than  to  re-en^hasize  the  many-electron 
nature  of  the  general  problem.  Because  of  the  reduced  dimensionality,  Coulomb  interactions 
among  ic-electrons  result  in  marked  electron  correlation  effects  which  render  independent 
particle  theories  mcomplete  for  description  of  nonlinear  optical  properties.  This  has  already 
been  demonstrated  in  the  comparison  of  experimental  and  theoretical  detenninations  of 
^jk(-2(0;(0,<i)),  where  calculations  involving  both  single  and  double-excited  configuration 
interactions  (SDCI)  result  in  good  agreement  with  experiment,  whereas  calculations  with  only 
single-excited  configuration  interactions  (SCI)  differ  by  as  much  as  50%.  In  the  case  of  Yiju 
for  linear  polyenes  to  be  discussed  below,  independent  particle  models  fair;  even  worse, 
resulting  not  only  in  different  magnitudes  but  also  opposite  sign  than  SDCI  calculations  and 
existing  experimental  data.  Because  independent  particle  results  are  not  only  quantitatively 
incorrect  but  also  miss  some  important  qualitative  features  of  and  Yijkl.  is  essential  to 
take  account  of  electron  correlations  through  multiple-excited  configuration  interactions  in 
descriptions  of  nonlinear  optical  processes. 
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A. 


One  Dimension;  G)njugated  Linear  Chains 


We  briefly  review  here  the  origin  of  the  microscopic  third  order  susceptibility 

in  prototype  conjugated  liiwar  chains  as  revealed  by  the  above  theoretical 
method. In  particular,  we  consider  the  linear  polyenes  which  are  hydrocarbon  chains 
in  which  each  carbon  site  is  covalendy  bonded  to  a  hydrogen  and  its  two  nearest  neighbor 
carbons.  The  remaining  valence  electron  of  each  carbon  atom  contributes  to  a  delocalized, 
strongly  onrelated  ic-election  distribution  altxtg  the  carbon  chain.  The  ground  stale  of  this 
system  is  a  spin-singlet,  broken-symmetry  state  in  which  the  carbon  lattice  possesses  a. single- 
bond/double-bond  altemadoa  We  have  performed  calculations  on  polyenes  ranging  in  length 
from  four  to  sixteen  carbon  sites  (N  «  4  to  16)  in  both  the  all-rron^  and  cis-transoid 
conformations. 

Although  the  self-consistent  field  calculation  of  the  grourxi  state  includes  all  of  the 
valence  shell  elections  for  each  atom  in  the  molecule,  Q  theory  needs  only  to  consider  rt- 
electron  orbitals  since  the  low-lying  exciutions  are  x  x*  transitions  and  for  conjugated 
systems  the  x-electron  contributions  to  ‘yiju(*<04;(i)i,oi>2,ct>3)  dominate  those  from  o-electrons. 
Because  the  polyenes  are  members  of  the  C2h  symmetry  group,  ail  of  the  x-electron  states 
must  possess  either  Ag  or  Bu  symmetry.  Within  the  ccntrosymmetric  C2h  group,  these  two 
symmetries  are  of  opposite  parity  leading  to  optical  dipole  selection  rules.  The  ground  state  is 
always  ^Ag,  and  therefore  the  'Bu  states  are  one-photon  allowed  transitions  observable  in  ihe 
linear  absorption  spectrum.  The  excited  ^Ag  states,  on  the  other  hand,  are  one-photon 
forbidden  but  two- photon  allowed.  Experimental  and  theoretical  studiesll^l '  of  one- 
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photon  and  two-photon  resonant  processes  in  finite  polyenes  have  shown  that  below  the  first 
optically  allowed,  dominant  1  ^Bu  state  is  located  a  strongly  electron  correlated  two-photon 
2iAg  state. 

Because  the  7C-electrt>n  distribution  is  easily  deformable  along  the  conjugation  axis  by  a 
perturbing  optical  electric  field  but  is  strictly  confined  in  the  two  transverse  directions,  it  is  to 
be  expected  that  the  dominant  tensor  component  of  Yijkl  will  be  Yxxxx  for  which  all  of  the 
electric  fields  are  aloig  the  axis  of  conjugan(».  We  define  the  x-axis  as  along  the  chain  and 

the  y-axis  as  the  perpendicular  direction  within  the  molecular  plane.  In  all  of  our  calculations 

{ 

of  varying  chain  length  and  ccmfnmation  of  linear  polyenes,  the  Ytxxx  conoponent  was  indeed 
found  to  be  far  larger  than  any  of  the  others.  As  an  example,  for  the  specific  case  of  trans- 
octatetraene  (rrons-OT)  (Rgure  la)  consisting  of  eight  carbon  sites  (N>8),  the  independent, 
nonvanishing  tensor  components  of  die  third  harmonic  susceptibility  Yjkl(*3<i>;Q),(D,0))  at  a 
nonresonant  fundamental  photon  energy  of  0.65  eV  907  ^m)  are  Yxxxx'^lS.S,  Ycyyx  = 

0.6,  Yyxxy^-^.  and  yyyyy^2  x  10*3^  esu.  All  components  involving  the  z-component 
necessarily  vanish  because  of  the  symmetry  constraints  of  the  re-electron  orbitals. 

The  symmetries,  energies,  and  relevant  transidon  dipole  moments  of  the  ten  lowest 
calculated  excited  states  of  mvts-OT  are  listed  in  Table  I.  All  153  of  the  calculated  states  are 
included  in  the  calculation  of  7ij*dt  i>ut  Table  I  saves  to  illustrate  the  important  points.  The 
columns  and  refer  to  the  x-components  of  the  transidon  dipole  moment  of  each 
state  with  the  ground  state  and  with  the  1  ^By  state,  respecdvely.  The  opdcal  selection  rules  are 
observed  in  the  vanishing  transidon  moments  g  for  all  the  ^Ag  states  and  M^,ib  ^ 

^Bu  states.  It  is  also  seen  that  the  l^Bu  state  has  by  far  the  largest  and  the  b^Ag  state  has 
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the  largest  As  a  result,  these  two  are  the  primary  contributors  to  y^uuu  as  will  be 
described  below. 

It  is  instructive  to  consider  the  individual  terms  in  the  sum  over  states  perturbation 
expansion.  Based  on  £q.  (S)  and  the  symmetry  selection  rules  described  above,  it  is  evident 
that  the  rc-electron  states  in  a  third  order  process  must  be  connected  in  the  series  g  — » ^Bu 
^  Ag  ^Bu  g.  For  centrosymmetric  structures,  third  order  processes  necessarily  involve 
virtual  transidtxis  to  both  one-photon  and  two-photon  states.  For  rronr-OT,  there  are  (153)^ 
terms  involved  in  the  summations  of  Eq.  (S).  However,  two  of  these  terms  are  an  order  of 
magnitude  larger  than  all  the  others  and  ctmstitute  70%  of  Yxxu-  The  remaining  terms  to  a 
large  extent  cancel  one  another  resulting  in  a  much  smaller  net  contribution.  In  both  of  the 
dominant  terms,  the  only  ^By  state  involved  is  the  dominant  low-lying  one-photon  l^By  x- 
electron  excited  state.  In  addition  to  its  low  excitation  energy,  this  state  is  important  because 
its  7.8  D  transition  dipole  moment  with  the  ground  state  is  more  than  three  times  larger  than 
any  other  ground  state  transition  moment  One  of  the  two  major  terms  comes  from  the  double 
sum  of  Eq.  (5)  with  both  of  the  intermediate  states  being  the  1  ^Bu.  In  the  case  of  the  double 
sum,  the  middle  intermediate  state  is  always  the  ground  state.  This  term  makes  a  negative 
contribution  to  Yuxx  below  resonance,  since  both  the  numerator  and  denominator  are  positive 
but  the  double  sum  has  an  overall  negative  contribution.  The  other  major  term  is  from  the 
triple  sum  witii  the  6^Ag  state  as  the  middle  intermediate.  This  state,  calculated  at  7.2  eV,  has 
a  large  transition  moment  with  1  ^By  of  13.2  D.  This  term  makes  a  positive  contribution  to 
Yxxxx  and  is  larger  than  the  first  leading  to  an  overall  positive  value  for  Yxxxx-  Importantly,  the 
6*  Ag  consists  of  60%  double-excited  configurations  indicating  it  is  highly  correlated.  SCI 


calculations  obtain  a  negative  value  for  Yxxju  because  they  do  not  adequately  describe  this  sute 
and  therefore  omit  its  large  contribution. 

Figure  2  displays  the  calculated  dispersion  curve  for  Yxxxx('3(D;CD,a),ci))  of  trans-Ol  as 
a  function  of  the  input  photon  energy.  The  first  resonance,  located  at  1.47  eV  (X=0.84  p.m) 
and  indicated  by  the  vertical  dash  in  the  figure,  is  due  to  the  3q)  resonance  of  the  1  ^Bu  state. 
The  second  singularity,  located  at  2.08  eV  (X=0.6  )im),  is  from  the  2co  rescmance  of  the  2^  Ag 
state.  As  seen  in  Eq.  (5),  the  ^Bu  states  v«ill  have  both  3(o  and  (O  resonances  in  third  harmonic 
generation,  whereas  the  ^  Ag  states  will  have  only  2ci}  resonances.  Other  third  order  processes 
have  different  resonance  selection  rules.  For  example,  in  DCSHG  both  die  ^Bu  and  ^  Ag  states 
(or,  more  generally,  both  one-photon  and  two  photon  states)  have  2a)  and  a>  resonances  while 
there  are  no  3a)  resonances.  It  should  also  be  noted  that  in  real  systems  natural  broadening  of 
electronic  states  wiU  prevent  divergence  at  the  resonances.  This  is  accounted  for  theoretically 
by  making  die  excitation  energies  in  Eq.  (5)  complex  with  the  inclusion  of  an  imaginary 
damping  term.  Unfortunately,  damping  constants  cannot  be  reliably  calculated,  and  one  must 
therefore  use  empirical  values. 

In  order  to  determine  the  influence  of  structural  geometry  on  nonlinear  optical 
prq}erties,  we  have  also  calculated  ‘Kju  for  the  cis-transoid  (cis)  conformation  of  polyenes 
(Figure  lb)  over  the  same  range  chain  lengths.  The  results  for  the  cis  conformations  are  in 
direct  analogy  to  those  for  tmns  indicating  that  geometry  is  not  a  predominant  factor.  For  the 
range  of  chain  lengths  we  have  considered,  we  find  that  the  transition  energies  of  the  1  ^  By  and 
2lAg  states  of  the  cis  confonnation  are  just  slighdy  red-shifted  from  the  values  for  irons  by 
0.02  to  0.10  eV.  For  example,  for  trtwu-OT  the  1  *  Bu  state  is  calculated  at  4.42  eV  and  the 
2^Agat4.16eV  while  for  cij-OT  the  values  are  4  37  eV  and  4.11  eV,  respectively. 


Furthermore,  just  as  for  trans,  the  dominant  tensor  component  of  Yijkl(-o>4;o)i, (02,0)3)  is  Yxxxx, 
and  the  most  significant  virtual  transitions  and  the  dispersion  are  essentially  the  same  as 
discussed  above.  The  mechanism  for  7ijki(-(04;(0i, (02,(03)  is  still  a  symmetry-dicuted  virtual 
excitation  process  involving  strongly  correlated  7t-electron  states.  The  calculated  dispersion  of 
Yxxxx(-3(o;(o,(0,co)  for  cij-OT  is  compared  to  that  of  trans-Ol  in  Figure  2.  The  two  curves  are 
very  similar,  and  the  most  significant  difference  is  the  smaller  nonresonant  value  in  the  case  of 
cis-OT.  Because  of  the  nearly  equal  excitation  energies,  the  first  two  resonances  occur  at  very 
near  the  same  energies  for  the  two  different  structural  conformations. 

The  similarity  between  the  two  conformaticms  is  further  emphasized  in  the  transition 
density  matrix  contour  diagrams  that  graphically  illustrate  the  electron  redistribution  upon 
virtual  excitation.  The  transition  density  matrix  pnn'  is  defined  through  the  expression 

<Hnn  >  *  -e  /  r  pnn'  (r)  dr  (6) 


with 


Pnn'  (ri) »  J  Vn  (ri,r2,...rM)  Vn  (ri.r2,.-rM)  dr2...drM  (7) 

where  M  is  the  number  of  valence  electrons  included  in  the  molecular  wavefunction  y. 
Contour  diagrams  for  Pnn’  of  the  ground  and  6*  Ag  states  with  the  1  *Bu  state  for  both  the  cis 
and  trans  forms  of  OT  are  compared  in  Figure  3,  where  solid  and  dashed  lines  correspond  to 
increased  and  decreased  charge  density.  The  contour  cut  is  taken  0.4  A  above  the  molecular 
plane  since  71  orbitals  vanish  on  the  atoms.  The  antisymmetry  of  Pnn'  is  a  requirement  for  a 
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one-photon  transition,  and  for  a  large  transition  dipole  moment  |inn'  the  contours  should  have 
well  separated  increased  and  decreased  densities  with  pardcularly  large  magnitudes  towards 
the  ends.  Of  course,  the  actual  charge  density  of  any  state,  given  by  Pnn.  must  be  symmetric 
because  of  the  centrosymmetry  of  the  molecule  so  that  all  states  will  have  vanishing  dipole 
moments.  In  Figure  3,  the  (a)  cis  and  (b)  rrans  virtual  g  1  ^Bu  transition  results  in  a 
somewhat  modulated  redistribution  of  charge  with  transition  moment  x-components  of  7.9  and 
7.8  D,  respectively.  For  the  1  ^Bu  ->  6^  Ag  transition,  however,  there  results  a  highly 
separated  charge  distribution  in  the  (c)  cis  and  (d)  trans  conformations.  Tlic  corresponding 
transition  mtxnents  are  12.0  and  13.2  D  for  the  cis  and  trans  cases,  respectively. 

There  is  very  little  experimental  data  available  for  Yiju(-Q>4;CDi, (02,0)3)  of  finite 
polyenes.  Our  calculations  cotxq)are  quite  well  with  the  experimental  values  that  have  been 
reported^^^^  for  two  of  the  shortest  polyenes.  Accounting  for  the  proper  fundamental  photon 
energy,  isotropic  orientational  averaging,  and  the  small  estimated  o-electron  contribution,  our 
results  are  in  good  agreement^  with  the  gas  phase  DCSHG  measurements  of  butadiene 
(N=4)  arxl  hexatriene  (N^ti)  by  Ward  and  Elliott  Comparison  of  calculated  results  with 
experimental  measurement  of  both  the  sign  and  magititude  of  y  is  obviously  an  important  test 
of  the  validity  of  a  theoretical  method,  but  it  has  all  too  often  been  neglected  in  recent  literature 

Analysis  of  the  chain  length  dependence  of  Yxxxx  for  the  trans  and  cis  conformations 
shows  that  the  two  are  further  unified  by  a  common  power  law  dependence  of  Yxxxx  on  L, 
defined  as  the  distance  in  the  x-direction  between  the  two  end  carbon  sites.  The  power  law 
exponent  was  found  to  be  4.6  ±  0.2  with  several  common  features  leading  to  this  rapid  growth 
of  Yxxxx  with  chain  length.  First,  the  excitation  energies  decrease  with  increased  chain  length 
Second,  the  magnitudes  of  transition  dipole  moments  along  the  chain  axis  increase  steadily 
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with  chain  length.  Thiid,  while  for  OT  and  the  shoter  chains  the  ncMilinear  susceptibility  is 
almost  entirely  composed  of  the  contributions  from  only  a  few  states,  longer  chains  have 
significant  conizibutions  from  an  increasingly  larger  number  of  both  ^Bu  and  ^  Ag  states. 

B .  Two  Dimensions:  Conjugated  Cyclic  Chains 

In  this  section,  we  extend  the  microscopic  description  of  Yijkl  for  one  dimension  to  two 
dimensions  and  consider  a  specific  example  within  the  class  of  conjugated  cyclic  structures 
known  as  annulenes.  As  a  major  case,  we  coitsider  the  planar  structure  of  cyclo-octatetraene 
(COT),  the  cyclic  analog  to  OT  with  N=8,  illustrated  schematically  in  Figure  4.  Although  the 
geometrically-relaxed  ground  state  of  COT  is  known  to  be  a  non-planar,  bent  structure,  we 
will  examine  only  the  planar  structure  of  COT  in  this  section.  The  purpose  is  to  explore  the 
effect  of  dimensionality  on  the  microscopic  feanires  of  Yjkl(*<)i>4;coi,(D2,ci>3).  Because  the 
inclusion  of  the  geometrically-relaxed  distortions  of  the  physically  observed  COT  structure 
unnecessarily  complicates  comparison  of  cyclic  structure  results  with  those  described  above 
for  linear  chains,  we  consider  only  die  planar  structure.  Analogies  and  contrasts  between 
linear  and  cyclic  structures  are  clearest  when  the  only  distinction  is  the  increase  in 
dimensionality  from  one  to  two  so  that  one  need  not  be  concerned  with  decoupling 
dimensionality  from  non-planarity  effects. 

For  linear  chains,  since  the  x-con^nent  of  the  transition  dipole  moments  is  much 
larger  than  the  y-  and  z-components,  the  isotropically  averaged  susceptibility  Yg  is  effectively 
determined  by  the  Yxxj«  component  while  all  components  involving  transverse  fields  make 
negligible  contributions.  This  reduces  Eq.  (4)  to 
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Yg  •  j  Yxxw  . 


(8) 


But  for  cyclic  structures,  since  the  in-plane  x-  and  y-directions  are  equivalent,  Yxx«  Yyyyy 
should  be  of  equal  magnitude.  Furthermore,  components  such  as  Jxxyy  ^  significant 

since  they  involve  these  two  directions  as  well.  The  corresponding  non-negligible  terms  of  Eq. 
(4)  are 


Yg  *  y  [  Yxxxx  ■*”  Yyyyy  ^  (Y«yy  "hyxy  +  fxyyx  +  +  Tiyxyx  +  T^nixy)]  •  (9) 


It  seems,  therefore,  that  one  ought  be  able  to  enhance  Yg  by  moving  from  linear  to  cyclic 
conjugated  structures  and  opening  pathways  for  new  con^onents  of  the  Yijkl  tensor  to 
contribute.  We  shall  demonstrate  below,  however,  a  nx>st  striking,  opposite  finding.  We 
show  that,  because  of  the  relevant  length  scales  involved  in  the  two  problems,  the  conjugated 
cyclic  structure  will  necessarily  have  a  smaller  yg  than  the  corresponding  ctxijugated  linear 
chain  with  an  equal  number  of  carbon  sites. 

Our  planar  model  COT  is  a  member  of  the  dihedral  Dah  symmetry  group  which  is 
non- Abelian  and,  hence,  has  two-dimensional  irreducible  representations  denoted  as  E  classes 
The  allowed  state  symmetries  for  it-electron  excitations  are  Aig,  A2g,  B ig,  B2g,  and  Eq.  Of 
these,  only  states  of  symmetry  are  one-photon  allowed  excitations  from  the  ^Aig  ground 
state.  The  ^Eu  states  are  doubly  degenerate  with  the  two  representations  related  by  a  7i/2 
rotation  about  the  z-axis  perpendicular  to  the  molecular  plane.  All  of  the  remaining  symmctne' 
listed  above  describe  nondegenerate,  two-photon  states.  A  typical  feature  of  conjugated  cyclic 
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molecules*  including  phthalocyanines  and  poqjhyrins,  is  the  existence  of  a  relatively  low 
frequency  absorption  in  the  visible  or  near  ultraviolet  and  a  higherirequency  absorption  deeper 
in  the  ultraviolet^^  This  feature  appears  in  our  model  COT  with  the  weak  low  frequency 
1  ^Eu  state  at  4.4  eV  and  the  much  stronger  higher  frequency  2^Eu  state  at  6.5  eV.  In  the  case 
of  phthalocyanines,  however,  it  is  always  the  low  frequency  band  which  is  stronger. 

In  Table  n  we  list  the  symmetries,  energies,  and  relevant  transition  dipole  mcnnents  for 
the  eight  lowest  calculated  excited  states  of  COT.  Again,  a  total  of  153  states  are  calculated. 
The  third  and  fourth  columns  list  the  x  and  y-components,  respectively,  of  the  transititm 
moment  between  a  given  state  and  the  ground  state  while  the  column  labeled  lin,2E 

component  of  the  transition  moments  of  that  state  with  the  two  degenerate  representations  of 
the  2^Eu  state;  the  y-components  are  given  in  the  column.  Although  the  x-  and  y- 

direcdons  are  equivalent  in  COT,  it  is  seen  in  the  Table  that  these  components  of  the  transition 
naoments  are  not  always  equal,  and  this  is  a  direct  result  of  the  double  degeneracy  of  the  Eu 
states.  Transition  naoments  involving  any  degenerate  pair  of  ^Ep  states  with  a  two-photon  state 
are  clearly  related  by  a  tc/2  rotation  or  x-»y,  y-»  -x.  Thus,  the  appearance  of  negative  signs  in 
some  of  the  transition  naoments  merely  reflects  the  choice  of  basis  in  Hilben  space  and  has  no 
absolute  physical  meaning.  By  choosing  an  appropriate  basis  for  the  degenerate  pair,  the 
magnimdes  of  the  x  and  y  components  can  be  made  equal,  although  this  is  not  necessary  as 
illusoated  in  Table  n. 

The  first  one-photon  doubly  degenerate  1  'Eu  state  has  very  small  transition  moments 
with  the  ground  state  and  a  correspondingly  small  oscillator  strength.  The  2^Eu  doubly 
degenerate  state,  on  the  other  hand,  has  a  maximum  transition  moment  of  4.6  D  in  this  x-y 
basis  and  is  a  large  oscillate  strength  one-photon  state  in  analogy  to  the  1  state  of  OT. 


I 


Thus,  the  transition  moments  of  the  two-photon  states  with  the  2^Eu  state  deteimine  their 
importance  for  Yijki(*0)4;(Oi,Ci)2,ci>3)  and  are  also  listed  in  Table  II  for  the  lowest-lying  states. 

Analysis  of  the  various  terms  in  the  summations  of  Eq.  (5)  for  COT  reveals  much 
similarity  to  the  linear  chain  analog.  As  in  OT,  there  is  one  dominant  term  in  the  double 
sununation  which  is  due  to  the  state  with  the  largest  oscillator  strength,  in  this  case,  the  2iEu 
state.  The  dominant  contributions  from  the  triple  sum  also  all  involve  the  2^Eu  state,  but  rather 
than  just  one  term,  there  are  several  significant  contributions  for  COT.  Each  one  is  smaller 
than  the  dominant  negative  term  from  the  2^Eu  state,  but  together  they  again  lead  to  posidve 
values  for  the  nonresonant  tensor  components  YjU-  The  significant  contribution  from  several 
terms  involving  different  intermediate  two-photon  states  is  amilar  to  what  was  found  for  the 
linear  polyenes  longer  than  OT. 

In  Figure  5,  the  transition  density  matrix  Pm'  is  shown  for  the  ground  state  with  the 
two  representations  of  die  2iEu  state.  The  ic/2  rotational  relationship  of  the  two 
representations  of  2iEu  is  also  clear  here  as  in  Table  II.  The  charge  redistribution  for  this 
transition  is  fairly  modulated  as  it  is  for  the  linear  chains.  In  the  particular  basis  shown  for  the 
representation  of  the  douUy  degenerate  2iEu  state,  the  transition  moments  alcmg  the  x  and  y 
directions  are  comparable  as  seen  in  the  Figure. 

The  calculated  dispersion  of  the  isotropic  third  order  susceptibility  Yg(-3(o;o),(0,CD)  is 
given  in  Figure  6.  The  noniesonant  value  is  quite  small  compared  to  the  linear  chains.  Very 
sharp  resonances  occur  at  1.47  eV  due  to  a  3o)  resonance  to  1  ^Eu  and  1.60  eV  due  to  a  2co 
resonance  with  2lAig.  Since  both  of  these  states  have  very  small  transition  moments  in  Table 
n,  they  only  make  significant  contributions  directly  on  resonance.  Calculations  including 
imaginary  damping  terms  to  account  for  the  finite  width  of  the  electronic  excitations  show  that 
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these  resonances  becozoe  completely  washed  out  and  the  dispersion  remains  flat  in  this  pan  of 
the  spectrum.  Experimental  dispersion  measurements,  therefore,  likely  would  not  observe 
significant  resonant  enhancement  until  the  3ci)  resonance  to  the  2^Eu  state  which  occurs  at  2. 16 
eV. 

The  calculated  values  of  Yiju(-3a);co,ca,C0)  for  COT  at  the  ncMiresonant  fundamental 
photon  energy  of  0.65  eV  are  Yxxxx=^*75  and  Y«yy=0.21  x  10*3^  esu.  Because  of  the  Dah 
symmetry.  Yxxxx^’l^yyy  and  Yuyy  ^  equal  to  all  other  components  that  involve  two  x- 
component  fields  and  two  y-component  fields.  This  equivalence  between  the  x  and  y* 
directions  is  in  contrast  to  the  linear  polyene  case  where  the  Yixxx  component  dominates  all 
others.  Furthermore,  the  Yxxxx  component  of  OT  (15.5  x  10*^  esu)  is  far  larger  than  any  of 
the  components  of  COT.  However,  since  COT  has  significant  components  in  both  the  x  and  y 
directions,  it  is  more  reasonable  to  conq)are  values  of  the  isotropically  averaged  susceptibility. 
For  COT,  Yg  *  0.38  x  10*36  esu  as  con^pared  to  3.4  x  10*36  esu  for  trans-OT.  As  in  the  case 
of  comparing  values  for  the  trans  and  cis  conformations  of  linear  polyenes,  here  also  we  must 
consider  the  actual  lengdi  scale  involved  in  the  problem. 

The  important  length  for  nonlinear  optical  responses  is  the  largest  length  over  which 
charge  can  be  separated  due  to  the  presence  of  an  optical  electric  field.  For  the  linear  polyenes, 
that  length  is  the  distance  along  the  conjugation  axis  separating  the  end  carbons.  For  COT,  the 
equivalent  length  is  between  two  points  on  either  end  of  a  diameter  of  the  ring.  Additionally, 
since  n-electron  motion  is  constrained  along  the  carbon  lattice,  the  relevant  length  is  one-half 

of  the  circumfence  of  the  molecular  structure.  This  distance  is  5.3  A  for  COT.  For 
comparison,  for  trans-KT  (N*6)  Yg  =  0.75  x  10' 36  esu  at  htD=0.65  eV  and  the  end-to-end 

length  L  is  6  A.  Thus,  for  these  two  very  different  structures  where  the  relevant  lengths  are 
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nearly  equal,  we  find  that  the  calculated  values  of  yg  are  also  comparable.  In  the  case  of  cyclic 
structures,  the  length  which  determines  the  magnitude  of  Yg  appears  to  be  one-half  of  the 
circumference,  and  this  is  a  general  Hnding  in  calculations  we  have  performed  for  other  size 
rings  which  will  be  published  separately. 

3.  RESONANT  MEASUREMENTS  IN  QUASI-TWO  DIMENSIONAL  RING 
I  STRUCTURES:  SATURABLE  ABSORPTION  AND  OPTICAL  BISTABILITY 

We  have  shown  in  the  previous  section  that  increased  dimensionality  finom  linear  to 
cyclic  chains  leads  to  a  decrease  in  the  isotropic  susceptibility  Yg  rather  than  an  enhancement 
as  one  might  first  expect  This  is  due  to  an  effective  reduction  in  the  length  available  for  tiie  x- 
electrons  to  respond  to  an  applied  optical  electric  field.  There  are,  however,  other  compelling 
reasons  to  investigate  conjugated  cyclic  structures  for  their  nonlinear  optical  properties. 
Primary  among  these  is  the  presence  of  optically  intense  low  and  high  fiequency  bands  in  the 
visible  and  ultraviolet  generally  observed  in  large  diverse  classes  of  related  ring  structures  that, 
in  addition,  exhibit  inqKxtant  secondary  material  properties  such  as  thermal  and  chemical 
stability  and  ease  of  fabrication  and  processing.  Among  such  classes  are  the  well-known,  large 
ring  porphyrin  and  phtiudocyanine  structures  which  exhibit  well-defined  intense  Q  and  Soret 
bands  in  the  visible  and  near  uv,  respectively.  Free  base  pcxphyrin,  for  example,  has  long 
been  identified  as  an  analog  of  18-armulene  which  possesses  nine  double  bonds.^^^^  Although, 
in  general,  we  expect  the  nonresonant  Xj^  of  ring  structures  not  to  be  very  large,  porphyrin- 
j  like  structures  provide  many  attractive  features  for  studying  resonant  Xj^  processes. 
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Optical  bistability  in  saturable  absorbers  has  long  been  under  intense  study, first, 
to  understand  bistable  phenomena  and  their  relationships  to  the  mediating  intensity  de¬ 
pendent  re&active  index  n2,  and,  second,  to  develop  thin  film  nonlinear  optical  devices. 
Based  on  the  microscopic  descriptions  discussed  in  the  previous  section,  we  have  selected 
phthalocyanine  structures  in  order  to  conduct  saturable  absorption  and  optical  bista¬ 
bility  studies  in  quasi-two  dimensional  conjugated  ring  structures.  We  have  succeeded  in 
designing  ultrathin  naphthalocysmine  films  possessing  large  resonant  nj  values  and  in  im¬ 
plementing  wide  aperture  Fabry-Perot  thin  film  etalons  for  optical  bistability  st  dies.  In 
addition  to  primary  nonlinear  optical  properties,  the  design  of  the  thin  film  etalons  incor¬ 
porated  important  secondary  material  properties:  (i)  large  saturable  absorption  with  zero 
unsaturable  background;  (ii)  saturable  absorption  centered  in  the  spectral  range  of  oper¬ 
ating  laser  diodes;  (iii)  thin  film  formation;  (iv)  spin  coatable  films  and  high  through-put 
fabrication;  and  (v)  photolithographic  properties  for  micron  size  pattern  features. 

A  first  major  example  is  the  silicon  naphthalocyanine  oligomer(SINC)  with  the  molec 
ular  structure  shown  in  Figure  7.  Both  pure  SING  and  SING  solid  solution  ultrathin 
were  directly  fabricated  by  spin  coating  techniques.  The  films  were  homogeneous  and  ; 
high  optical  quality  with  typical  thickness  of  80nm  and  420nm  for  pure  dye  and  solid  S'  I  ; 
tion,  respectively.  Figure  8  shows  the  linear  absorption  spectra  of  the  pure  SING  filmi  5.  . 
line)  and  a  solid  solution  film  of  SING  in  PMM.\(dashed  line).  The  optical  spectrum  • 
hibits  an  intense  Q-band  centered  in  the  near  infrared  at  810nm  (1.53  eV)  for  the 
dye  and  774nm  (1.60  eV)  for  the  solid  solution,  with  a  linewidth(FWHM)  of  0.1  eV  a., 
linear  absorption  coefficient  a  of  1x10*  cm''  The  molecular  symmetry  of  SING  is  I ' ; 
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and  the  doubly  degenerate  Q*band  is  ansilogous  to  the  l^Ai^  — » ir-electron  transition 
of  COT  described  in  the  previous  section. 

We  first  present  a  physical  model  for  the  optical  excitations  of  the  thin  film  and  then 
present  saturable  absorption  studies.  A  Maxwell-Bloch  two-level  model  is  used  for  the 
dependence  of  the  saturation  threshold  power  on  pulse  duration,  and  the  experimental 
results  support  this  simple  picture.  For  laser  frequencies  near  the  absorption  maximum, 
we  find  an  effective  na  of  1x10'^  cm^/kW.  Based  on  the  saturable  absorption  behavior, 
electronic  absorptive  optical  bistability  at  nanosecond  time  scales  is  then  observed  in  thin 
film  Fabry-Perot  etalons.  Bistability  behavior  is  studied  as  a  function  of  longer  pulse 
duration  purposely  to  distinguish  thermally-induced  bistable  effects. 

A.  Glassy  Polymer  Dye  Model  and  Saturable  Absorption 

In  this  section  we  present  a  physical  model  for  optical  excitations  in  the  SINC  ul- 
trathin  films.  The  films  obtained  by  spin  coating  do  not  possess  either  positional  or 
orientational  long  range  order.  Because  of  this  random  distribution  of  SINC  molecules  in 
the  films,  the  excitations  are  strictly  on-site  7r-electron  transitions,  which  have  a  typical 
absorption  coefficient  a  of  10*  cm~'.  The  on-site  7r-electron  optical  excitations  in  an  iso¬ 
lated  molecule  have  an  intrinsic,  temperature  independent  natural  linewidth  on  the  ord^r 
of  0.1  to  1  GHz  with  corresponding  radiative  decay  lifetime  of  1  to  10  ns. 

The  large  linewidth  of  the  Q-band  observed  in  Figure  8  is  primarily  due  to  inhonv. 
geneous  broadening.  This  site  broadening  in  the  thin  film  has  a  Gaussian  shape  which 
a  consequence  of  the  statistical  distribution  of  resonance  frequencies  of  the  optical  cent**'- 
due  to  a  variation  in  local  environment  in  the  polvmer  matrix.  Within  the  inhomogene.. 
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bro&dened  Gftussian  envelope  of  the  Q-band  are  a  series  of  narrow  homogeneous  broadened 
resonances  where  the  characteristic  temperature  dependence  of  the  linewidth  depends  on 
the  microscopic  broadening  mechanism.  The  linewidth  of  the  homogeneous  broadening 
in  amorphous  media  can  be  accounted  for  using  the  TLS(two  level  system)-glass 

models  introduced  by  Anderson,  Halperin  and  Varma  which  have  been  fairly  successful 
in  accounting  for  many  of  the  physical  properties  (e.g.  specific  heat,  thermal  conductivity, 
ultrasonic  absorption)  of  disordered  systems  such  as  glasses  and  polymer-like  matrices. 
The  Hamiltonian  of  the  glassy  polymer  dye  thin  film  on  a  substrate  can  be  written  as  ^34] 


H  =  Hii  +  ^  Ha  (lOfl) 

Ha  =  co'J'o^o  +  +  \e<t*  +  Yi  Wn  +  \) 
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a  a 
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where  Ho  is  the  noninteracting  Hsuniltonian  of  the  molecule  (the  difference  ci  —  co  corrr 
spending  to  the  Q-band  absorption),  the  TLS,  and  the  phonons;  H12  is  the  electrostaur 
dipole  interaction  between  the  molecule  and  the  TLS;  and  H23  is  the  interaction  beiwrrr 
the  TLS  and  the  phonons,  i.e.,  the  strain  field  <  is  coupled  to  the  TLS. 

Each  homogeneous  line  under  the  Gaussian  envelope  is  approximated  by  a  Loren t::.i 
function,  and  the  width  is  related  to  the  temperature  dependent  population  and  phase  ■ 
laxation  rate  of  the  excited  state  by  Fourier  transformation.  The  amorphous  mater  < 
surrounding  an  optical  site  can  be  considered  as  an  ensemble  of  noninteracting  TLS 


flopping  between  two  eigenstates  as  they  emit  or  absorb  acoustic  phonons.  The  homoge¬ 
neous  line  broadening  in  amorphous  media  comes  from  the  dipole  coupling  of  the  optical 
sites  with  the  TLS’s,  and  the  line  width  can  be  expressed  in  terms  of  TLS  the  lifetime  (or 
flip-flopping  rate),  the  temperature  dependence  of  which  is  determined  firom  the  coupling 
between  the  TLS  and  the  strain  field  manifested  as  acoustic  phonons.  In  this  way,  the 
line  width  of  optical  sites  in  the  amorphous  media  is  predicted  to  have  T^*^^  (0<  S  <1) 
dependence,  and  the  absolute  magnitude  can  be  obtained  once  all  the  physical  parameters 
for  the  TLS  are  known.  For  various  organic  molecules  in  different  polymeric  matrices,  the 
theoretically  predicted  temperature  dependence  of  the  homogeneous  linewidth  is  in  good 
agreement  with  spectral  hole  burning  data.  According  to  photochemical  hole  burn¬ 

ing  data  on  porphyrins  and  phthalocyanines  in  various  glasses  and  polymeric  matrices, 
the  homogeneous  linewidth  of  the  naphthalocyanine  oligomer  at  room  temperature  is  es¬ 
timated  to  be  between  10  and  100  GHz.  For  SING  films,  comparison  of  the  homogeneous 
linewidth  (100  GHz)  with  the  inhomogeneous  width  (0.1  eV)  shows  that  the  Gaussian 
envelope  contains  on  the  order  of  ten  thousand  Lorentzian  broadened  resonances, 

Saturable  absorption  behavior  within  the  Q-band  of  the  thin  pure-SINC  and  SIN( 
solid-solution  filmii  was  investigated  as  a  function  of  pulse  duration  and  wavelength  ot  ih>- 
incident  beam.  A  methane  Raman  cell  pumped  by  fundamental  or  doubled  Nd:YAG  out  pu  ■ 
was  used  to  produce  Stokes  radiation  near  the  peak  maximum  of  the  Q-band;  at  813nm  '  ' 
a  pure  SING  thin  film  and  770nm  for  a  SINC  solid  solution  thin  film.  The  methane  Raiu  > 
cell  pumped  by  a  Quantel  Nd:YAG  picosecond  laser  system  provided  30GHz  bandwidth 
ps  pulses  of  200  y.i  per  pulse  that  were  easily  focused  to  give  light  intensities  exceedn  v 
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MW/cm^.  Nanosecond  pulses  of  30  GHz  bandwidth  with  an  energy  of  20  /xJ  per  pulse  were 
obtained  in  the  same  way  by  using  a  Quanta-Ray  DCR  Nd:YAG  laser  system  to  provide 
an  incident  intensity  of  500  kW/cm^.  Figure  9  shows  the  incident  intensity  dependence  of 
the  absorption  oL  for  30  ps  (circles:  pure-SINC,  triangles:  SING  solid'solution)  and  10  ns 
(squares:  pure-SINC)  pulses,  respectively.  The  change  in  the  absorption  was  reproducible 
through  many  cycles  of  increased  and  decreased  incident  light  intensity.  Data  points  for  10 
ns  are  more  scattered  owing  simply  to  power  fluctuations  in  the  Nd:YAG  pump  laser.  In 
the  case  of  the  pure  SING  film,  if  the  light  intensity  was  increased  more  than  the  maximum 
of  the  data  points,  i.e.,  30  MW/cm^  for  30  ps  pulses  and  200  KW/cm^  for  10  ns  pulses,  an 
irreversible,  simple  melt  phase  change  was  observed  in  the  film.  But  for  the  solid  solution 
film,  no  irreversible  change  was  observed  up  to  an  incident  intensity  of  4GW/cm^,  and 
furthermore,  the  saturation  behavior  was  almost  identical  to  that  of  the  pure  SING  film. 

The  solid  lines  in  Figxire  9  are  least  squares  curves  to  a  Bloch-type  saturable  absorp¬ 
tion 

where  I©  is  the  threshold  power  for  the  saturation,  ooL  gives  the  low  intensity  linear 
absorption,  and  obL  is  the  unsaturable  background  absorption.  Importantly,  asL 
found  to  be  zero  for  both  the  pure  SING  and  solid  solution  films.  In  the  case  of  the  p '  r-- 
SINC  film,  the  threshold  powers  for  saturation  were  lOOMW/cm^  and  440kW/cm^  for 
and  10ns  pulses,  respectively.  The  threshold  power  for  10ns  pulse  is  lower  than  thai  ■ 
30ps  pulse,  but  this  trend  saturates  at  or  near  the  10/xs  scale.  Importantly,  for  the  Sl'«> 
solid  solution  thin  film,  the  threshold  power  for  10 ps  pulse  was  also  lOOMW/cm^,  the  '  < 
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value  as  that  of  the  pure  SINC  film.  The  fact  that  the  saturable  absorption  behavior  is 
identical  for  both  pure  SINC  and  solid  solution  film  is  quite  striking,  and  can  be  accoimted 
for  by  the  absence  of  any  positional  and  orientational  order  between  molecular  sites  in  the 
thin  filtn  phase;  thus,  there  is  no  phase  coherence  between  optical  sites.  From  this  we  cam 
conclude  that  the  on-site  ir-electron  excitations  of  the  Q-band  in  individual  molecular  sites 
are  responsible  for  the  large  resonant  nonlinear  optical  response. 

The  pulse  duration  dependence  of  the  saturation  absorption  threshold  power  is  well 
explained  in  terms  of  a  two-level  optical  system  interacting  with  the  resonant  incident 
light^^°^(^^l .  The  dynamic  transmission  can  be  described  by  two  coupled  rate  equations, 

=  -/(i,l)»ln(i,()  -  {N,  -  n(*,l)}l  -  (12) 

cn  T, 

=  -I{x,t)(r[{No  -n(i,0}  -n(x,<)]  (13) 

where  No  is  the  total  number  of  molecules  per  unit  volume,  n{x,t)  is  the  number  of 
molecules  in  the  excited  state  at  depth  z  amd  time  t,  a  is  the  absorption  cross  section  per 
molecule,  and  t,  is  the  radiative  decay  time  of  the  excited  state.  Now  let  the  incident  pulse 
shape  f{i)  be  normalized  by  the  incident  light  intensity, 

/(o,<)  =  /o/(o 

As  the  light  pulse  paisses  through  the  absorbing  molecule,  the  pulse  is  distorted  as  wet! 
attenuated.  The  change  in  the  pulse  shape  can  be  expressed  in  terms  of  a  transrrus- 
function  T(x,t)  at  depth  z  and  time  t,  defined  as 

/(z, t)  =  /o/i  n D X.  () 
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The  two  equations  (12)  and  (13)  can  be  reduced  to  one  equation  for  the  transmission 

T{x,t): 


dliiT(x,t) 

dii/r,) 


+  In  T{x,t)  =  2<rr,/o/(f)(l  —  T{x,t))  +  In  To 


(16) 


where  /q  is  the  incident  light  intensity,  To  is  the  low  intensity  transmission,  and  f{t) 
is  the  incident  light  pulse  profile.  From  numerical  simulation  of  Eq.(16),  the  saturation 
threshold  intensity  I,  can  be  obtained  for  a  given  pulse  width  and  excited  state  lifetime 
T,.  Furthermore,  the  ratio  of  saturation  threshold  for  two  given  pulse  widths  can  be  used 
to  uniquely  determine  r«.  Using  230  as  the  ratio  of  threshold  powers  at  30  ps  and  10 
ns  as  determined  from  our  experiment,  we  find  the  excited  state  lifetime  to  be  5  ns,  in 
good  agreement  with  previous  independent  relaxation  measurements  on  related  structures 

[30l[42l_ 


From  the  saturable  absorption  measurement,  we  can  also  estimate  the  nonlinear 
refractive  index  n2.  The  nonlinear  optical  susceptibility  x(‘*')  for  a  Bloch  type  system  is 
given  by 


=  _ )  ,i: 

^  ^  +//// 

where  oo  is  the  linear  absorption  coefficient,  A  is  the  detuning  (u;  -  u;o)/r,  and  I,  , 
the  saturation  power.  The  nonlinear  refractive  index  n2  is  the  derivative  of  the  intent  i  . 
dependent  refractive  index  n  with  respect  to  the  inensity  I,  so  for  this  model 
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From  this  expression,  we  obtain  a  value  of  nj  for  the  pure-SINC  thin  film  of  1  x  10  ^cm 
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To  fully  ch&racterize  the  nonlinear  absorption  properties  of  the  pure  SINC  film,  we 
measured  the  dispersion  of  the  Q-band  nonlinear  absorption.  A  Spectra-Physics  380B  cw 
ring  dye  laser  using  LDS  821  dye  and  tunable  between  790  nm  and  860  nm  was  used 
as  the  light  source  with  typical  peak  powers  of  150  mW.  The  cw  dye  output  of  5  MHz 
bandwidth  was  chopped  by  an  acousto-optic  modulator  (AOM),  so  that  10  /xs  pulses  were 
produced  at  a  low  duty  cycle  (10~'*)  to  eliminate  any  heating  effects.  A  light  intensity 
of  approximately  80  kW/cm^  was  obtained  by  focusing  the  laser  output  directly  onto  the 
sample  with  a  microscope  objective  lens. 

The  dispersion  of  the  nonlinear  absorption  of  the  thin  film  is  compared  in  Figure  10 
with  the  linear  absorption  spectrum.  It  is  clear  that  the  change  in  absorption  is  maximum 
on  resonance  and  decreases  in  moving  away  from  the  peak.  The  dashed  lines  are  a  least 
squares  fit  to  a  Gaussian  envelope  model  for  the  linear  absorption,  and  the  dispersion 
was  derived  from  the  Kramers-Kronig  relations.  The  dotted  lines  are  the  fitted  curves 
for  the  nonlinear  absorption  and  dispersion.  We  can  sec  that  the  maudmum  change  of 
refractive  index  as  obtained  from  Eq.(18)  occurs  at  776  nm,  and  the  maximum  n2  at  776  nm 
results  from  a  linear  superposition  of  the  refractive  index  changes  of  each  homogeneous! \ 
broadened  resonamce  within  the  inhomogeneously  broadened  envelope. 

B.  Optical  Bistability 

SINC  thin  films  can  be  employed  as  the  nonlinear  medium  in  Fabry-Perot  etalons  '  • 
optical  bistability  studies.  The  etalon  was  formed  by  spin  coating  a  wide  area  (2 
in  diameter)  thin  film  of  approximately  80  nm  thickness  on  the  front  mirror  of  two  hic-  • 
reflecting  dielectric  mirrors  of  a  Burleigh  RC-110  Fabry-Perot  interferometer.  The  • 


■pectral  range  of  the  cavity  was  adjustable  and  set  at  12.5  GHz.  Initial  cavity  detuning 
was  adjusted  by  varying  the  high  voltage  applied  to  a  PZT  annular  ring  that  held  the 
output  mirror  (the  PZT  had  total  motion  of  2.21  /im/1000  V),  and  when  this  output 
mirror  was  scanned,  a  Fabry-Perot  interference  pattern  with  a  finesse  of  2  was  obtained 
with  the  sample  in  the  cavity. 

The  optical  bistability  experimental  layout  (Figure  11)  was  basically  the  same  as 
the  arrangement  used  for  measuring  the  dispersion  of  saturable  absorption.  The  laser 
output  was  partially  reflected  by  a  pellicle  beam  splitter  to  a  reference  silicon  detector, 
and  the  transmitted  beam  was  tightly  focused  on  the  thin  film  which  had  been  spin  coated 
directly  onto  the  front  mirror.  The  output  fringes  were  collimated  and  passed  through  a 
500  nm  pinhole,  and  a  sample  silicon  detector  monitored  the  intensity  at  the  central  pan 
of  the  Bull’s  eye  interference  fringe.  Silicon  PIN  photodiodes  were  used  for  reference  anu 
sample  arms,  and  the  output  intensity  versus  input  intensity  was  directly  displayed  on  a  ; 
GHz  oscilloscope  operating  in  the  x-y  mode  and  externally  triggered.  At  nanosecond  timr 
scales,  initial  measurements  were  attempted  at  813  nm  using  the  output  from  the  Ndi’t'A^  < 
pumped  methane  Raman  cell  sis  the  laser  source,  but  later,  more  stable  amd  smooth  pul'P' 
at  799  nm  from  a  single-mode  pulsed  TiiAljOj  laser  pumped  by  a  frequency  dour  • 
Nd:YAG  laser  were  used.  Figure  12  shows  the  bistable  hysteresis  behavior  observed  a  ■ 

40  ns  pulses  with  intensity  160  kW/cm^  at  799  nm  were  incident  on  the  nonlinear  F  P 
found  that  for  zero  cavity  detuning  (cp- O.Oi  the  bistable  behavior  W2is  the  largest,  a 
when  the  cavity  was  detuned  (<p=7r),  the  shape  of  the  hysteresis  curve  changed  and 
effect  became  much  smaller.  In  each  case,  the  data  were  reproducible  through  many  <  , 


of  increased  and  decreased  incident  light  intensity.  Further,  the  data  were  idso  reproducible 
when  focussing  at  different  areas  of  each  fUm  sample.  This  kind  of  initial  cavity  detuning 
dependence  is  typical  of  absorptive  OB,  which  results  from  saturable  absorption  and  a 
subsequent  change  in  the  loss  of  the  F-P  cavity,  enabling  higher  transmission  than  is 
possible  at  low  light  levels.  However,  the  greatly  reduced  effects  observed  for  nonzero 
cavity  detuning  indicate  a  very  small  dispersive  contribution. 

The  incident  laser  wavelength  was  moved  to  780  nm  corresponding  to  near  optimum 
conditions  for  possible  dispersive  contributions  to  the  bistability.  The  same  measurements 
were  repeated,  and  under  all  conditions  no  bistable  behavior  was  observed  (Figure  13).  The 
absence  of  any  bistable  behavior  at  780nm,  where  a  large  dispersive  effect  is  expected  from 
the  Kramers-Kronig  relation  analysis,  means  simply  that  the  relatively  broad  Q-band  is 
not  one  single  homogeneously  broadened  line  but  an  inhomogeneously  broadened  envelope 
for  many  homogeneous  lines,  and  supports  the  microscopic  picture  of  the  glsissy  polymer 
SING  film.  Thus,  near  maximum  saturable  absorption  at  810  nm,  naphthalocyanine  thin 
films  exhibit  primarily  absorptive  optical  bistability  at  nanosecond  time  scales.  The  above 
analysis  is  based  on  steady  state  conditions  for  optical  bistability  because  of  the  relevant 
time  scales  (pulse  width:  40ns;  relaxation  time:  5ns;  and  cavity  roundtrip  time:  In^ 
but  transient  effects  may  also  play  some  role  in  the  observed  hysteresis  behavior,  whi- 
continues  to  be  examined  in  on-going  studies. 

Since  purely  absorptive  optical  bistabiUty  is  rarely  observed,  a  systematic  sene-^ 
studies  of  possible  thermally-induced  refractive  index  changes  in  the  thin  films  were  pr- 
formed  as  a  function  of  increased  pulse  duration  from  10'*  to  2  seconds.  The  laser  sour  • 


used  in  these  studies  was  the  AOM  chopped,  Spectra-Physics  cw  ring  dye  laser  described 
above.  Figure  14  shows  an  example  of  the  optical  bistability  observed;  the  pulse  width  was 
2.0  sec  and  the  duty  cycle  20%.  When  the  cavity  was  tuned  on  resonance  with  the  incident 
light  wavelength  (<^=0.0,  lower  right  in  Figure  14),  no  hysteresis  loop  was  observed,  indi¬ 
cating  the  lack  of  an  absorptive  bistability  effect,  and  this  was  true  even  when  the  dye  laser 
output  frequency  was  resonant  with  the  linear  absorption  peak  of  the  sample.  However,  for 
a  negative  ^  (-1.0,  left  in  Figure  14),  bistable  hysteresis  with  a  coimterclockwise  circulation 
was  observed,  and  for  a  positive  ^  (0.75,  upper  right  in  Figure  14),  the  bistable  hysteresis 
had  clockwise  circulation.  With  the  pulse  width  and  the  resonant  wavelength  near  the 
linear  absorption  peak  of  810  nm,  the  bistable  behavior  is  due  to  a  thermally-induced 
dispersive  intensity  dependent  refractive  index  change.  Most  importantly,  this  behavior  is 
distinctly  different  from  the  fast  pulse  bistability  observed  at  nanosecond  timescales. 

For  both  electronic  absorptive  and  thermal  dispersive  optical  bistability,  the  hysteresis 
loops  do  not  show  a  complete  switching  on  and  switching  off.  The  reason  for  this  is  tha( 
eventhough  the  resonant  nj  value  of  the  SINC  thin  film  is  furly  large,  the  film  thickness 
(80  nm)  is  so  small  that  the  nonlinear  phase  shift  experienced  by  the  light  inside  the  etalon 
is  not  large  enough  for  complete  switching  to  occur.  Thicker  SINC  containing  films  ha\<* 
been  prepared  in  our  continuing  studies. 

4.  CONCLUSION 

Electron  correlation  effects  markedly  determine  the  virtual  electronic  excitation 
cesses  and  nonlinear  optical  properties  of  the  low  dimensional  structures  of  conjukia  ■ 
linear  and  cyclic  chains.  Microscopic  understanding  of  the  the  origin  of  the  large. 
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fast  nonresonaat  Xijfc(“‘^s;u;i,wa)  and  xl*fcj(~‘*>4;u>i,wa,u>s)  of  these  materials  is  obtained 
through  many-electron  calculations  of  the  molecular  susceptibilities  and 

wijWjjWj)  and  direct  comparison  with  available  experimental  results.  Details  of 
the  mechanism  for  7»i*i(~‘*>4;wj,u>2,44;s)  of  one  dimensional  conjugated  linear  chains  have 
been  explained  in  terms  of  virtual  transitions  among  highly  correlated  ir-electron  states, 
and  contour  diagrams  of  transition  density  matrices  pnn'  provided  direct  illustration  of 
the  most  significant  virtual  transitions.  The  dependences  of  7iasi(~ur4;a/i,a>],u;3)  on  chain 
length  and  structural  conformation  have  been  examined,  and  conformation  was  found  to 
be  only  significant  in  as  much  as  it  affects  the  actual  physical  length  of  the  chain.  Results 
for  the  all- irons  and  cis~transoid  conformations  of  polyenes  are  unified  by  a  power  law 
dependence  of  the  dominant  tensor  component  7„.«(~3w;w,u>,u>)  on  the  chain  length  L 
with  an  exponent  of  4.6i:0.2. 

A  theoretical  study  of  7i;*i(— u;4;u>i,u>j,u>3)  in  a  two  dimensional  conjugated  cycU' 
structure  allowed  examination  of  the  role  of  dimensionality  in  the  nonlinear  optical  prop¬ 
erties  of  conjugated  systems.  The  origin  of  7ijij(— u>4;u>j,u;3,u>3)  in  this  case  was  found  to 
be  very  similar  to  that  for  the  one  dimensional  case,  but  the  magnitudes  of  the  various 
components  of  jijki  are  all  smaller  than  the  component  of  the  one  dimensional  chaiu 

with  an  equal  number  of  carbons.  More  important,  the  isotropically  averaged  third  ord'*r 
susceptibility  7^  is  also  much  less  for  the  two  dimensional  structure  as  compared  to  the  'i  • 
dimensional  chain.  Our  analysis  shows  that  this  is  a  result  of  a  reduction  in  the  effect 
length  over  which  ir-electrons  can  respond  to  an  applied  optical  electric  field.  For  t  ■ 
dimensionad  cyclic  structures,  this  length  is  one-half  of  the  circumference  of  the  nnc 
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compared  to  the  full  end-to-end  chain  length  in  one  dimensional  structures.  As  a  general 
rule,  two  dimensional  conjugated  cyclic  structures  are  therefore  expected  to  have  smaller 
nonresonant  values  than  their  one  dimensional  analogs. 

The  presence  of  large  oscillator  strength  x-electron  bands  in  the  visible  and  near 
ultraviolet  characteristic  of  conjugated  quasi-two  dimensional  structures  provides  attrac¬ 
tive  conditions  for  studying  resonant  processes,  especially  in  phthalocyanine-related 
structures  which  also  possess  outstanding  secondary  material  properties.  We  have  demon¬ 
strated  that  such  two  dimensional  structures  can  be  designed  as  wide  area,  spin  coatable 
nonlinear  optical  films  and  that  these  can  be  represented  as  a  microscopic  composite  sys¬ 
tem  of  molecular  optical  sites  in  a  TLS-glass  random  medium.  The  resonant  nonlinear 
optical  properties  of  homogeneously  broadened  lines  contained  in  the  inhomogeneously 
broadened  Gaussian  envelope  were  directly  investigated  by  standard  saturable  absorption 
measurements.  These  results  showed  for  the  SINC  films  that  saturable  absorption  occurs 
at  the  peak  maximum  of  the  low  frequency  Q-  band  near  810  nm  with  essentially  zero  un- 
saturable  absorption  background  and  that  on-site  rr-electron  excitations  of  the  Q-  band  m 
individual  molecular  sites  are  responsible  for  the  large  resonant  nonlinear  optical  response 
n2  of  1x10"^  cm*/kW  at  810  nm. 

A  cueful  study  of  optical  bistabihty  in  thin  film  etalons  demonstrated  that  near 
the  peak  maximum  of  the  saturable  Q-band  centered  at  810  nm,  SINC  thin  films  exhih;: 
primarily  electronic  absorptive  optical  bistabihty  at  fast  time  scades  (nanoseconds  i 
long  timescales  (seconds),  thermally-induced  dispersive  bistability  occurs  which  is  eniir'*  . 


different  and  distinct  from  the  observed  fast  time  beliavior. 


The  present  results  represent  a  first  report  on  studies  centering  on  conjugated  quasi- 
two  dimensional  structures.  Father  theoretical  and  experimental  studies  of  larger  ring 
structures  and  micron  thick  films  of  naphthalocyanines  for  their  resonant  and  nonresonant 
behavior  have  been  completed  and  will  be  reported  separately. 
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Table  Q.  The  symmetiies,  energies,  and  selected  transition  dipole  naoments  of  the  calculated 
low-lying  states  of  cyclo-ocatetraene. 


Symmetry 

Energy  (eV) 

4n^(D) 

I'Aj, 

2.24 

0.00 

0.00 

-1.60,  2.40 

-2.40,  -1.60 

2'Ai, 

3.19 

0.00 

0.00 

0.14,  0.09 

-0.09,  0.14 

I'E, 

4.41 

0.05 

-0.03 

0.00,  0.00 

0.00,  0.00 

i'e. 

4.41 

0.03 

0.05 

0.00,  0.00 

0.00,  0.00 

I'Bj, 

5.21 

0.00 

0.00 

-0.33,  0.50 

0.50,  0.33 

I'Bxi 

5.94 

0.00 

0.00 

1.03,  0.68 

0.68,  -1.03 

2>E„ 

6.48 

4.58 

-3.04 

0.00,  0.00 

0.00,  0.00 

2‘e, 

6.48 

3.04 

4.58 

0.00,  0.00 

0.00,  0.00 
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Figure  1.  Schematic  diagrams  of  the  molecular  structures  for  (a)  ail-trans  and  (b)  cis-transoid 
polyenes. 

Figure  2.  Calculated  dispersion  of  Yxxxx(*3o);Q},ci},a})  for  trans-OI  (solid  curve).  The  vertical 
dashes  locate  the  3<o  resonance  to  the  l^Bu  state  and  the  2(o  resonance  to  the  state.  The 
dotted  curve  shows  the  analogous  caknilared  dispersion  for  cis-OT. 

Figure  3.  Transition  density  matrix  diagrams  for  (a)  cis-OT  and  (b)  trans-Orr  ground  states 
with  the  l^Bu  state  and  (c)  cis-OT  and  (d)  trons-OT  6^Ag  states  with  the  l^Bn  state.  In  both 
cases,  diese  are  the  most  inqxxtant  virtual  transitions  for  liitxxx(‘<A4;t>oi>tb2*c»3)-  The 
corresponding  x-components  of  the  transiticm  dipole  moments  are  7.9, 7.8, 12.0,  and  13.2  D, 
respectively. 

Figure  4.  Schematic  diagram  of  the  molecular  structure  of  planar  cyclo-octatetraene  (COT). 

Bgure  5.  Transitioo  density  matrix  diagrams  for  the  ground  state  of  COT  with  each  of  the 
two  rqnesentations  of  the  doubly  degenerate  2^Ey  sute.  The  it/2  rotational  relationship 
between  the  two  representations  is  clear.  In  this  basis  for  the  2iEu  degenerate  pair,  the  x  and 
y-components  of  the  transition  dipole  moments  are  4.6  and  -3.0  D  for  (a)  and  3.0  and  ^  6  D 
for  (b),  respectively.  Other  choices  for  the  2*Eu  basis  will  result  in  different  values  for  the 
individual  oon^nents  of  the  transidon  dipole  moments,  but  the  magnitude  of  the  transidon 
dipole  moment  vector  must  remain  invariant 
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Figure  6.  Calculated  dispersion  of  the  isotropically  averaged  third  order  susceptibility 
79(— 3w;a>,u;,u;)  for  COT. 

Figure  7.  Schenaatic  diagram  of  the  moleculair  structure  of  silicon-naphthalocyanine 
oligomer  (SINC). 

Figure  8.  Linear  absorption  spectra  of  pure  SINC  (solid  curve)  and  SINC  solid  solution 
(dashed  curve)  thin  films.  The  large  oscillator  strength  Q-band  absorption  peaks  at  1.53 
eV  (810  nm)  in  pure  SINC  film  and  1.60  eV  (774  nm)  in  SINC  solid  solution  film. 

Figure  9.  Saturable  absorption  data  and  least  squares  fit  to  Eq.  (11).  Squares  are  data 
for  pure  SINC  film  at  pulse  width  rp=  10  ns;  circles,  pure  SINC  film  at  30  ps;  and 
triangles,  SINC  solid  solution  film  at  rp=:  30  ps.  aL  is  reflection  calibrated  and  the  low 
intensity  limit  is  aLs  1.40. 

Figure  10.  Dispersion  of  saturable  absorption  between  1.45  eV  and  1.55  eV  measured  with 
10  /iS  pulse.  The  solid  line  is  the  linear  absorption  curve.  The  dashed  line  is  the  Gaussian 
function  least  squares  fit  curve  of  the  linear  absorption,  and  the  dotted  line  is  a  similar 
least  squares  fit  of  the  nonlinear  absorption  The  curves  with  maxima  occtirring  at  A=  77^ 
nm  correspond  to  the  real  parts  of  the  refractive  indices  obtained  from  the  Krame.'s-Kronie 
relations. 

Figure  11.  Schematic  diagram  of  optical  bistability  experiment  for  the  pulse  widths  10'*  ' 
2.0  sec  (AOM  :  acousto-optic  modulator,  PH  :  pinhole,  PD  :  photodiode).  For  nanosecor  i 
ptilse  width  mesisurements,  a  single-mode  Ti:AljOj  laser  replaced  the  Ar"*"  laser,  ring 
laser,  and  AOM  configuration. 


Figure  12.  Absorptive  optical  bistability  observed  for  polymer  dye  etalon  with  pulse  width 
Tp=  40  ns  and  wavelength  A=799  nm.  In  (a)  the  initial  cavity  detuning  is  set  at  0,  and 
in  (b)  the  initial  cavity  detuning  4>o  is  set  at  ir. 

Figure  13.  No  optical  bistability  observed  for  polymer  dye  etalon  with  pulse  width  Tp=  40 
ns  and  wavelength  A=780  nm  for  initial  cavity  detuning  set  at  either  0  or  ir. 

Figure  14.  Thermal  dispersive  optical  bistability  observed  for  polymer  dye  etalon  with 
pulse  width  rp=  2  sec  and  wavelength  A=810  nm.  The  lower  right  curve  is  for  the  initia 
cavity  detuning  0.0,  the  left  curve  for  ^o=  -1*0,  and  the  upper  right  c\irve  for  ^o= 
0.75. 
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II.  FREQUENCY  DEPENDENCE  OF  THE  OPTICAL  KERR  EFFECT  AND  THIRD 
ORDER  ELECTRONIC  NONLINEAR  OPTICAL  PROCESSES  OF  ORGANIC 
SYSTEMS 

I.  INTRODUCTION 

The  optical  Kerr  effect^ (OKE)  is  a  four  wave  mixing  process  in  which  an  intense 
optical  pump  beam  of  frequency  (Oj  induces  depolarization  in  a  probe  beam  of  frequency  (O2 
because  of  the  third  order  nonlinear  optical  polarization  P.^^(ci)2).  In  an  isotropic  medium, 
pfk)  is  expressed  as 

Pp^ojj)  “  Ej(C0,)Ek((0,)*Ei(t02),  (1) 

where  XjjJ  is  the  third  order  nonlinear  optical  susceptibility  and  E(cD|)  and  £(0)2)  are  the  pump 
and  probe  Helds,  respectively.  As  extensively  discussed  in  the  review  by  Hellwarth,^  at 
frequencies  well  above  molecular  vibraticm  and  rotation  modes  but  betow  electronic  resonances, 
has  contributions  due  to  both  electronic  excitations  and  lattice  displacements  which  in  the 
Bom-Oppenheimer  approximation  is  expressed  as  the  sum  of  the  quantum  averages  and 

Dijki* 

where  and  are  the  third  order  electronic  and  lattice  contributions,  respectively,  and 
is  a  local  field  factor  relating  micioscoptc  and  macroscopic  quanitities ;  both  tenns  are  assumed  to 
be  fast  relative  to  experimental  time  scales.^'^  In  such  a  frequency  regime,  because  of  higher 
frequency  electronic  resonances  nearby,  the  dominant  frequency  dependent  term  is  the  ultrafasi 
electronic  contribution,  and,  hence,  frequency  dependent  measurements  of  the  total  OKE 
response  can  allow  direct  determination  of  the  resonant  electronic  and  nonresonant  contnbunor  \ 
to  X^^\  This  is  similar  to  the  approach  of  Levenson,  et  al.,^®*''  where  the  nonrcsonjr 
contribution  was  determined  by  comparing  OKE  measurements  with  complementary  Raman 
scattering  data  (Raman-induced  Kerr  effect  or  RDCES).  The  present  work  aims  to  determine  v 


separate  resonant  electronic  and  nonresonant  contributions  from  OKE  data  over  a  large  frequency 
range. 

The  purpose  of  this  paper  is  to  provide  an  experimental  and  theoretical  study  and  review 
that  enables  frequency  dependent  OKE  results  to  be  obtained  from  optimized  pump-probe 
techniques.  In  Section  2,  the  probe  beam  depolarization  is  theoretically  analyzed  in  terms  of 
experimentally  measured  quantities  with  special  attention  paid  to  several  standard 
approximations.  Section  3  describes  the  experimental  pump-probe  design  and  the  formalism 
required  for  the  data  analysis.  Section  4  presents  experimental  results  for  the  frequency 
dependent  measurements  of  X^j^(-(D2:ci>p-C0j.0)2)  liquids  benzene  and 

nitrobenzene.  In  Section  5,  electronic  and  lattice  mechanisms  are  discussed.  A  simple  model  is 
presented  to  describe  the  dispersion  of  sum-frequency  and  difference-frequency  resonances  in 
benzene  and  nitrobenzene  in  terms  of  their  linear  absorption  spectra  and  dipole  selection  rules, 
and  to  obtain  the  resonant  and  noniesonant  contributions  to  Xiju(-o)2;co^,-<Dj,o)2)  for  the  rwo 
liquids. 

2.  THEORETICAL  APPROXIMATIONS 

In  this  section,  explicit  account  is  given  for  approximations  and  assumptions  required  for 
optimizing  the  experimental  design  of  OKE  measurements  using  pump-probe  techniques.  The 
plane  wave  limit  is  considered  under  the  conditions  that  the  probe  beam  waist  is  small  lelaove  o 
the  pump  waist  and  that  the  effects  of  the  pump  on  the  probe  are  the  same  everywhere  on  ihe 
transverse  profile  of  the  probe.  Also,  both  the  probe  and  pump  beam  profiles  have  waists  whcs« 
diffraction  and  self-focusing  lengths  are  much  longer  than  the  sample  interaction  length,  insur.r  < 
that  the  puiz^  and  the  probe  beams  at  small  crossing  angles  will  result  in  a  uniform  effect  on  -  f 
probe  beam  throughout  this  interaction  region.  Since  the  peak  power  of  the  pump  beam  i.'i 
present  study  is  enough  to  produce  a  measurable  refractive  index  change,  unfocused  beams  i-r 
used  allowing  a  long  interaction  length  that  minimizes  extraneous  nonlinear  effects  such  as  :  ~  ^ 


from  the  sample  cell  windows.  The  pump-probe  conditions  are  reviewed  in  terms  of  Maxwell 
equations  in  the  following  discussion. 

In  an  isotropic  medium,  the  nonlinear  interaction  between  an  intense  spatially  delocalized 
pump  field,  E^,  and  a  weak  probe  field,  E^,  mediated  by  an  optical  Kerr-like  medium  with  a 
scalar  third  order  susceptibility,  is  expressed  by  the  nonlinear  Maxwell 

equation: 

-V2E  -  -  ^  ^  [  e^E  +  (E  E)El,  (3) 

where  the  total  electric  field  is  composed  of  tlw  pump  and  probe  fields, 

E  »  E|  +  Ej.  (4) 

With  terms  retained  only  to  first  order  in  E^  in  Equation  (3),  the  pump  intensity  acts  as  a  source 
term  in  the  wave  equation  for  the  probe. 

For  plane  waves,  Ej  and  Ej  are  given  by 

Ej »  j  [Aje*^  +  C.C.1  X,  (5) 

and 

E]  a  y  [Aje^^e^'^cosSx  +  sin0y)  +  c.c.J,  (6) 


where 
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-  k^z  -  ©it  +  0(z), 


(7) 


(8) 


0(z)  is  the  intensity  dependent  phase  shift  of  the  pump,  6  the  angle  that  the  polarization  vector  of 
the  probe  makes  with  that  of  the  pump  at  the  sample  entrance  (z^O),  5(z)  a  phase  shift  induced  in 
the  probe  by  the  pump,  and  and  A2  may  be  assumed  constant  if  the  beams  approximate  plane 
waves  in  the  measured  region.  Equation  (S)  assumes  that  ail  the  effects  of  the  probe  on  the  pump 
are  negligible  so  that  the  pump  renoains  linearly  polarized  at  all  points  in  the  interaction  region. 
Note  that  the  wide  pump  waist  implies  that  the  phase  shift,  5(z),  is  independent  of  the  transverse 
coordinate  resulting  in  a  uniform  phase  shift  in  the  transverse  plane.  Without  this  assumption, 
the  measurement  would  require  knowledge  the  probe  phase  everywhere  in  the  detection  plane, 
or  the  use  of  a  deconvolution  analysis,  that  would  greatly  complicate  data  interpretation.  If  <)Kz) 
or  8(z)  is  complex,  then  for  each  quantity  the  real  part  corresponds  to  an  intensity  dependent 
phase  shift  while  the  imaginary  part  describes  intensity  depeiulent  absorption.  In  either  case, 
there  results  a  ptdaiization  change  in  the  probe  beam. 

Equation  (3)  can  be  Fbuiier  analyzed  to  obtain  the  equations  satisfied  by  the  pump  and 
probe  beams.  We  assume  that  the  induced  nonlinearity  is  small  enough  that  the  phase  changes 
negligibly  on  wavelengtii  scale.  The  phase  change  of  the  probe  beam  due  to  the  pump  is  then. 


5(z)  - 

efcaj) 


The  output  to  input  ratio  of  the  probe  beam  in  a  Kerr  shutter  is 


•  sin^20)sin^  [ 
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where  L  is  the  interaction  length  in  the  sample.  In  the  limit  of  small  nonlineaiities  we  can  sou  r 
Equation  (10)  for 
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(U) 


It  is  significant  to  note  diat  Equation  (1 1)  is  only  valid  when  the  argument  in  brackets  in  Equation 
(10)  is  small  as  required  by  the  slowly  varying  phase  approximation;  otherwise,  the  formalism 
fails  and  higher  order  effects  will  also  contribute,  contrary  to  the  assurtq)tions  of  small  refractive 
index  change.  If  8(z)  is  complex  (8(z)  =  8^(z)  +  i8i(z)),  the  solutions  follow  a  parallel 
development  where  sin^(8/2)  is  replaced  with  sin^(8j^/2)+sinh^(8j(2),  which  in  the  small 
argument  limit  becomes  (8^/2)^  +  (8{/2)^.  Equadon  (11)  can  then  be  generalized  by  replacing 
with  -  ((XjV  + 

At  this  point,  it  is  worthwhile  to  discuss  the  approximations  we  have  employed  in 
arriving  at  Equation  (11)  and  how  well  the  experimental  parameters  obey  the  assumptions.  For  a 
Gaussian  pump  beam  of  initial  beam  waist  a,  V-E  is  negligible  in  the  nonlinear  wave  equation 
under  the  following  conditions:^^ 


X2 

Tt^a^f* 
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and 


We  further  require  the  probe  beam  to  be  contained  by  the  pump 
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(14) 


(15) 
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where  a,  is  the  probe  beam  waist  Here,  f  is  the  normalized  pump  beam  width  (f  »  w  is  the 
beam  width  as  a  function  of  z,  X  the  wavelength,  r  the  radial  distance  from  the  symmetry  axis  in 
a  plane  transverse  to  the  propagation  direction,  z  the  distance  along  the  propagation  direction,  and 
is  the  nonlinear  radius  of  curvature  of  the  beam  defined  as 
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In  the  experiments  described  here,  we  have  X  >  1  iim,  a  s  0.5  cm,  f~l,  z~200  cm,  lAjI^  »  3.4S  x 
10^  ergs/cm^,  esu,  and  =  1.5.  Thus,  typically  XV)i^V~4.1xl0'^  and 

z^/R^~0.05  satisfy  conditions  (12)  and  (14),  respectively. 

In  summary  for  the  plane  wave  analysis  to  be  valid,  the  sample  length  must  be  much 
shorter  than  both  the  diffraction  length  and  self-focusing  lengdi  scale,  and  the  probe  beam  must 
be  totally  contained  within  the  central  potion  of  the  pump  beam.  Moreover,  the  phase  shift  of 
the  probe  introduced  by  the  pump  must  be  small  on  the  scale  of  the  probe  beam  wavelength. 


3.  EXPERIMENTAL 

3.1  APPARATUS 

The  experimental  layout  is  shown  in  Figure  1.  A  oxxle-locked  Nd:YAG  laser  ((^aniel 
401C)  is  used  as  the  punq>  and  a  frequency  stabilized  argon  ion  laser  (Spectra  Physics  165)  or 
frequency  stabilized  ring  dye  laser  (Spectra  Physics  380B)  is  used  as  the  probe.  The  Nd:  Y  AG 
pump  laser  beam  is  split  into  four  pieces  by  beamsplitters  BSl  and  BS2,  so  that  two  of  the 
resulting  beams  are  used  to  pump  the  sample  (S)  and  reference  (R)  Kerr  cells  defined  bv 
polarizers  Pl/Al  and  P2/A2  respectively.  The  other  two  Nd:YAG  beams  are  optically  deta^ed 
with  stepping  motor  controlled  near-infrared  mirror  pairs  (YMl/MDMl,  YM2/MDM2)  and 
recombine  with  the  Kerr  signal  at  the  optical  gate  (KDP).  Similarly,  the  cw  argon  ion  or  dvr 
probe  laser  beam  is  split  between  the  sample  and  reference  arms.  A  weak  reflection  from  the  fuJ 
wave  plate  HWPl  illuminates  the  photodiode  that  gates  the  electronic  data  acquisition  sysiem 


The  probe  beam  frequency  stabilization  prevents  mode  hopping,  which  insures  that  the  beam  is 
continuous  wave  on  picosecond  time  scales.  High  power  near  infrared  interference  filters  are 
used  variously  as  beamsplitters  or  minors,  and,  when  possible  45-45-90  prisms  steer  the  beams 
by  total  internal  reflection.  Glan-Taylor  prisms  are  used  in  the  Kerr  cell  because  of  their  high 
rejection  rates,  and  an  interferometrically  aligned  quartz  sample  cell  prevents  probe  beam 
depolarization,  refraction  and  scattering.  The  optically  gated  signal  is  detected  with  a 
photomultiplier  tube  (RCA  C31034-A02),  and  the  subsequent  current  pulse  is  in  turn  sent  to  a 
discixxtinator,  a  CAMAC  gated  integrator  (LeCroy  2249W),  and  finally  to  a  mictocomputer  (DEC 
PDF  1 1/23)  for  analysis. 

Since  only  nonfocused  laser  beams  are  used,  requirement  (13)  is  satisfied.  Present 
experiments  used  NdtYAG  lasers  of  3-7  mm  width  and  cw  dye  and  Ar*^  lasers  of  1-2  mm  width 
consistent  with  Equation  (15).  Also,  the  central  part  of  the  measured  beam  profiles  are  consistent 
with  the  parabolic  approximation  of  self-focusing  calculations  except  for  the  small  dip  at  the 
origin  (Figure  2).  If  the  beam  crossing  angle  is  adjusted  between  two  runs,  this  dip  could  result 
in  a  ±8%  change  in  the  probed  pump  power  if  the  area  of  beam  overlap  changes.  When  relative 
values  of  n^  are  measured  between  two  samples,  this  can  result  in  a  total  uncertainty  in  this  rano 
of  about  16%.  As  long  as  no  adjustments  are  made  in  the  overlap  regions,  the  relative  frequency 
dependence  of  the  third  order  susceptibility  of  a  sample  will  be  correct  within  the  statisncai 
uncertainty  of  each  data  point.  Small  nonlinear  phases  (5  «  1)  were  obtained  by  employing  a  2 
cm  long  sample  cell  and  Nd:YAG  lasers  delivering  a  total  peak  power  of  Pj  >  1x10^^  erg/sec 
This  places  a  limit  on  the  lowest  measurable  value  of  n2,  which  in  the  present  study  is  10'^^  esu 

Any  Kerr  cell  experiment  must  address  the  problem  of  stray  light  getting  through  the 
crossed  polarizers  and  possibly  saturating  the  detector.  This  background  is  particularly 
significant  in  the  current  measurements  since  a  cw  probe  laser  is  employed  that  results  m  an 
average  signal  to  noise  ratio  of  10*^.  The  use  of  a  pulsed  pump  and  pulsed  probe  would  pcnr.n 
gating  via  synchronization  but  would  also  present  difficulty  in  separating  true  nonlinear  opncai 


processes  from  possible  pulse-noise  interactions  and  effects;  thus,  we  use  a  pulsed  pump/cw 
probe  arrangement  with  an  optical  gating  technique. 

A  KD*P  parametric  mixing  crystal  is  used  to  combine  the  Kerr  cell  signal  at  frequency  0)2 
and  a  portion  of  the  Nd:  YAG  laser  output  at  frequency  coj  to  create  ultraviolet  output  at  the  sum 
frequency  eu3sC0j+(i)2  by  noncollinear  phase  matching.  The  unwanted  light  at  cOj  is  filtered 
(spatially  and/or  spectrally)  allowing  only  the  ultraviolet  sum  frequency  light  to  reach  the 
detectors.  Only  when  the  pulse  of  light  emerging  from  the  Kerr  cell  is  synchronized  with  the 
pump  laser  arriving  at  the  KD*P  crystal  will  any  signal  be  detected  (temporal  filtering).  If 
operation  of  the  Kerr  cell  is  blocked  or  if  the  two  pulses  are  not  synchronized,  no  output  is 
obtained  from  the  optical  gate.  Accurate  synchronizadon  of  picosecond  pulses  is  achieved  using 
a  motorized  near-infrared  mirrcv  for  optical  delay.  Hence,  the  Kerr  cell  combined  with  the 
optical  gate  provides  a  highly  efficient  means  for  reliably  determining  whether  true  OKE  signal  is 
being  detected.  The  optical  gate  up-converts  10  to  20  percent  of  the  signal,  and  the  efficiency  of 
the  gate  is  monitored  as  a  function  of  the  input  parameters  such  as  laser  fm{uency  and  intensity 
crossing  angle  in  the  crystal,  etc.  When  properly  used,  the  optical  gate  is  an  essential  component 
in  experiments  using  cw  probes  by  providing  high  signal  to  noise  ratios  and  straightforward 
deconvolution  of  the  Kerr  signal  Furthermore,  the  optical  gate  can  be  easily  modified  for  nme 
resolved  studies  through  auto/cross  correlation  measurement  techniques.*^ 

3.2  DATA  ANALYSIS 

A  data  acquisition  system  employing  a  sample  and  reference  arm  measures  the  iirr.e 
integrated  power  of  each  laser  shot  during  the  electronic  gate.  The  integrated  charge  spectrum  ot 
at  least  1000  laser  shots  is  tracked  over  time  and  stored  and  the  centroid  of  the  dismbunon 
(average  number  of  photons  per  event)  is  calculated.  A  typical  pulse  height  distribution  of  si  gr.xi 
of  about  10,000  events  shows  the  expected  Poissonian  form  (Figure  3).  The  centroid  •. 
expressed  as  the  time  integrated  power  of  the  instantaneous  signal  power,  Pom(t),  and  tor  j 
single  pulse  the  contribution  is, 
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(17) 


C(ca2)  =  G(ca2)  / 


-At« 


where  0(0)2)  detection  gain  and  At  is  the  width  of  the  gate  in  the  detection  electronics.  For 
At  much  larger  than  the  temporal  pulse  width.  €(0)2)  can  be  expressed  as 


C(co2)  »  pL(“2)G(o»2)  /  q(0dt. 


(18) 


and  the  measured  pump  energy,  E. 
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where  P^(<«>2)  and  A®  are  the  peak  power  and  peak  field  values,  respectively,  A  the  area  of  the 
pump  beam,  c  the  speed  of  light,  and  q(t)  and  r(t)  axe  dimensionless  pulse  shape  functions 
having  their  respective  maTimtim  values  of  unity  at  M).  Substituting  Equations  (18)  and  ( W  > 
into  equation  (1 1),  we  obtain 
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A  reference  arm  is  used  to  eliminate  the  effects  of  pump  beam  power  fluctuations.  V- e 
reference  atm  can  use  either  the  same  probe  as  the  sample  arm,  or  a  separate  laser  source.  I n  ; ne 
former  case,  it  is  necessary  to  choose  a  reference  material  with  negligible  dispersion  in  rr 
the  probe  tuning  range.  The  latter  case  is  easier  to  realize  experimentally,  but  the  former  has 
advantage  of  eliminating  pump  and  probe  power  fluctuations  simultaneously.  In  the  dual 
configuration  with  all  variables  relevant  to  the  reference  arm  labeled  as  prime,  the  relative  4 
from  Equation  (20)  is  given  as 
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where  CCcOj)  and  C'(0),)  are  the  measured  centroids,  co^  the  reference  frequency,  nfcOj)  the  linear 
refractive  index  at  cOj,  and  Pin(0)2)  is  monitored  periodically.  In  deriving  Equation  (21),  it  is 
assumed  that  the  reference  arm  probe  frequency  is  fixed.  The  detection  gain,  G((02)  is  written  as 


*  «Gpmt(“2)Go,(“2)*  (22) 

where  Gg^IcOj)  is  the  gain  of  the  optical  gate,  Gp^fooj)  the  total  gain  of  the  photomultiplier  tube, 
and  the  constant  a  depends  on  the  experimental  geometry.  Since 
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where  MfcOj)  is  a  crystal  tuning  parameter  for  fQD*P,  Equation  (22)  yields 
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The  bracketed  term  in  (24)  is  constant  throughout  the  experiment  in  the  dual  probe 
configuration,  and  the  last  two  factOTS  can  be  measured  as  follows.  If  we  define  Pj^  as  the  Ken- 
effect  power  emerging  from  the  crossed  polarizers,  P,  the  light  passing  through  the  crossed 
polarizers  due  to  depolarization  from  sample  imperfections,  and  P^  the  power  permitted  to  pass 
through  by  slight  uncrossing  of  the  analyzer,  then  the  total  power,  P^,  reaching  the  optical  gaie 
is 


out 


p  +p 

s  u 


■•■p. 


I  I  3 


(25) 


The  purpose  of  uncrossing  the  polarizers  by  a  small  amount  is  to  tap  out  a  known  amount  of 
probe  light  in  the  same  spatial  direction  and  position  as  the  signal  beam.  Operation  of  the  optical 
gate  with  this  known  light  level  is  then  used  to  calibrate  the  gate  as  follows.  A  cw  power  meter 
measures  P,  when  the  polarizers  are  crossed  and  the  pump  beam  is  off;  +  P^  is  measured 
when  the  polarizers  are  slightly  uncrossed  and  the  pump  beam  is  off.  Similarly,  the  crossed  and 
uncrossed  Kerr  measurements  with  the  pump  on  yield  centroids  corresponding  to  light  at  the 
probe  frequency  entering  the  optical  gate  with  power  P,  +  Pj^and  P,  +  Pu  respectively. 
The  gain  of  the  optical  gate,  G(ai2),  is  the  ratio  of  the  measured  signal  centroid  corrected  for  the 
background  light  detected  by  the  PMT  and  the  known  input  intensity  and  is  given  by 
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where  $(0)2)  is  the  signal  centroid  &om  the  sample  normalized  to  the  reference  when  the 
polarizers  are  crossed,  and  $^{<02)  is  the  signal  centroid  when  the  polarizen  are  slightly 
uncrossed.  Substituting  Equations  (26)  and  (24)  into  Equation  (21),  we  finally  obtain 
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The  values  of  quantities  in  Equation  (27)  are  directly  experimentally  measurable. 

The  single  shot-to-shot  referencing  of  the  present  experiment  should  be  contrasted  with 
signal  averaging  techiuques.  Clearly,  in  either  case,  random  fluctuations  can  be  reduced  hv 
taking  moie  data;  however,  in  terms  of  systematic  errors,  the  methods  are  different  Many  opnes 
experiments  have  the  tendency  to  drift  over  time  and  can  be  quantified  by  a  time  constant,  t^.  .hit 
is  defined  as  the  elapsed  time  between  optimal  alignment  and  drift  to  the  e*^  point  of  the  measunro 
signal  relative  to  some  quasi-static  baseline.  In  the  single  shot-to-shot  experiment  the  refererKco 
signal  is  mostly  sensitive  to  inequalities  between  the  drift  rates  of  the  two  arms  of  the  expen  me 
(to  *tQ ).  Laser  pulse  shape  and  amplitude  fluctuanons  can  systematically  skew  the  data  by  v  ir-.x 
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of  their  unknown  distributions.  With  the  single  shot  referencing  method,  shape  and  amplitude 
changes  are  divided  out  for  each  individual  pulse,  thus  eliminating  any  need  for  knowledge  of  the 
shape  and  amplitude  distributions. 

Often,  in  time  domain  studies,  a  tempmal  pulse  shape  is  simply  assumed,^’^^*^®  Even  if 
temporal  profiles  are  known,  theoretical  calculations  sometimes  require  approximating  the  profile 
by  some  simple  function.  In  the  signal  averaging  case,  the  sample  is  measured  either  before  or 
after  the  reference,  hence  the  results  obtained  from  reversing  the  orders  can  be  averaged  to 
minimize  drift  error.  The  final  referenced  signal  in  the  signal  averaging  case,  however  will 
include  a  quotient  of  sums  of  pulse  amplitudes  and  shapes  from  the  two  separate  experimental 
runs.  The  value  of  this  quotient  cannot  be  predicted  without  some  a  priori  knowledge  of  the 
system  behavior,  even  though  the  general  behavior  of  the  system  can  be  very  well  quantified.  ^  ^ 
Hence,  the  single  shot  technique  minimizes  errors  owing  to  laser  fluctuations  induced  by  natural 
factors  such  as  temperature  change,  building  vibrations,  and  power  line  surges  as  well  as 
inherent  instabilities. 

4.  RESULTS 

4. 1  THERMAL  EFFECTS 

In  a  pump-probe  experiment  absorptive  beam  heating  effects  can  cause  changes  in  ihe 
sample  density  and  hence  in  the  sample  index  of  refraction.  Probe  heating  can  lead  to  thenr.ai 
defocusing  of  the  sigrud  that  results  in  spurious  variations  in  the  optical  gate  efficiency  as  well  as 
changes  in  the  pump  characteristics,  while  pump  heating  can  cause  probe  beam  defocusing  ar  j 
erroneous  apparent  contributions  to  the  total  measured  Kerr  coefficient  Although  heating  etTev  '.s 
are  small  at  fast  time  scales,  they  can  still  be  significant 

As  an  example,  the  thermal  response  time  of  nitrobenzene  was  measured  by  monitor  -  < 
the  amount  of  light  passing  through  the  sample  and  a  pinhole  centered  on  a  beam  of  waveler  £  ' 
X  a  488  nm  and  power  of  100  mW  as  a  function  of  time  after  exposure.  Typical  results 
given  in  a  log  plot  (Figure  4)  suggesting  a  response  nme  of  about  one  second.  A  scanr  -  < 
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pinhole  was  used  to  measure  the  spatial  beam  protile  as  a  function  of  beam  power  as  shown  in 
Figure  5.  The  difference  between  the  incident  beam  width  and  the  thermally  defocused  beam 
width  in  nitrobenzene  was  used  to  determine  the  thermal  contribution  to  the  Kerr  coefficient 
(Figure  6)  that  was  determined  from  self-defocusing  theory'*  to  be  about  0.1  esu. 

The  thermal  Kerr  coefficient,  n^,  depends  on  the  dme,  t,  after  initial  illumination  by  a  step 
function  source  and  the  relevant  material  response  time,  tg, 

nj  -  t^2(l-e-‘^),  (28) 

where  is  die  long  time  thermal  Kerr  coefficient  With  a  20ns  pump  pulse  at  wavelength  := 
488  nm,  the  thermal  Kerr  contribution  is  about  2x10'^  esu.  At  a  typical  pump  wavelength  of 
1064  nm,  the  absorption  coefficient  is  over  two  orders  of  magnitude  smaller  resulting  in  a 
thermal  contribution  of  less  than  2x10*"  esu.  This  upper  limit  is  nearly  the  electronic  and 
orientational  n2  of  single  ring  aromatic  molecules  like  nitrobenzene  and  the  orientational  n^  of 
molecules  like  CS2.  To  insure  that  thermal  effects  remain  negligible,  pump  pulse  durations  of 
less  than  l(X)ps  are  required.  Similarly,  the  long  time  thermal  effects  of  the  cw  probe  will  be 
negligible  for  probe  powers  less  than  a  few  mW.  Pump  heating  will  not  directly  cause  probe 
depolarization,  but  can  influence  its  magnitude  through  the  corresponding  refractive  index 
change. 

4.2  EFFECTS  OF  BEAM  PARAMETERS 

The  probe  power  leaving  the  Kerr  cell  is  given  by  Equation  (10).  The  signal  is  maximum 
when  the  angle  between  the  pump  and  probe  polarization  vectors  is  45°  and  the  measurej 
dependence  of  the  signal  (Figure  7)  on  the  polarization  angle  of  pump  beam  confirms  that 
signal  varies  as  sin^29).  The  pump  polarization  is  adjusted  to  maximize  the  signal  and  is  tl  x  (^d 
for  all  further  measurements. 


1 

In  the  limit  of  small  pump-induced  phase  shift,  the  power  leaving  the  Kerr  cell,  is 
related  to  the  probe  power,  Pj„,  and  the  pump  power,  Pj^p,  by 

j 

Po,.  -  (29) 

where  bis  a  power  independent  constant  calculable  from  equadon  (10).  We  note  that  expression 
I  (29)  is  strictly  valid  if  the  pump  induced  phase  change  is  small  and  if  the  response  time  of  the 

material  nonlinearity  is  fast  compared  with  the  laser  pulse  width;  this  is  always  true  for  the  non- 
resonant  electronic  contribudon.  If  the  response  time  of  the  material  is  comparable  to  the  laser 
pulse  width,  the  signal  will  be  a  convoludon  of  the  pump  laser  pulse  shape  and  the  response 
funcdon  of  the  mariinm.  Equadon  (23)  shows  that  the  gain  of  the  t^dcal  gate  is  propordonal  to 
the  product  of  and  Pp„^  since  both  the  gate  and  the  pump  beams  are  split  from  the  same 
source.  The  measured  signal  is  thus  propmdonal  to  the  cw  power,  Pj^,  and  to  the  cube  of  the 
pump  power,  Ppy^.  for  a  purely  third  order  process  and  small  pump-induced  phase  change. 

At  low  pump  powers,  the  signal  is  cubic  in  the  pump  power  as  expected  for  a  Kerr  effect. 
At  higher  pump  powers  (over  40  MW/cm^),  the  signal  is  no  longer  strictly  cubic  in  Ppu„,p  arid 
exhibits  behavitv  consistent  with  saturadon,  large  phase  shifts,  or  higher  mder  processes.  If  the 
pump  power  is  varied  only  in  the  Kcir  cell,  the  signal  should  be  quadratic  in  Pp„„^  and  this  was 
shown  for  pump  powers  less  than  40  MW/cm^.  The  dependence  of  the  signal  on  the  power  ot 
the  gate  pulse  in  the  KD*P  crystal  is  linear.  This  validates  the  calibration  of  the  optical  gate  as 
calculated  by  Equation  (23).  The  signal  is  also  linear  in  the  probe  power,  showing  th<r 
absence  of  thermal  defocusing  effects  up  to  10  mW  at  X  »  590  nm.  The  long  tail  of  the  broad 
absorption  ^ak  of  the  first  strong  excited  state  in  nitrobenzene  at  260  nm  explains  the  lower 
defocusing  threshold  of  about  1  mW  at  488  nm  as  opposed  to  the  higher  value  at  590  nm. 

The  experimental  beam  powers  are  set  within  limits  consistent  with  the  theorct;^  a. 
approximations  to  permit  the  signal  to  be  related  to  the  third  order  susceptibility.  In  the  prtx  e  ■> » 
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of  changing  pump  and/or  probe  frequencies,  the  thresholds  for  thermal  self-defocusing  and 
system  con^onent  saturation  must  be  monitored  and  adjusted  accordingly. 


4.3  DISPERSION  MEASUREMENTS 

Nitrobenzene  and  benzene  are  chosen  as  model  Kerr  liquids  because  the  wavelength  of 
their  ultraviolet  electronic  transitions  satisfy  conditions  for  sum>resonance  of  the  pump 
wavelength  at  1.06|zm  with  the  visible  cw  probe  wavelength.  The  pure  liquid  nitrobenzene 
absorption  spectrum  is  characterized  by  a  broad  strong  absorption  baiul  centered  at  4.6  eV  with 
weaker  bands  between  3  and  4  eV.  Gas  phase  absorption  spectra  show  that  there  should  be  a 
broad  weak  absorption  band  at  3.8eV  in  the  pure  liquid  appearing  as  a  small  shoulder  at  the  onset 
of  the  4.6  eV  band.  In  benzene,  the  first  strong  absorption  band  is  centered  at  4.86eV  and  a 
weak  band  in  the  near  ultraviolet  is  located  at  3.6S  eV.  The  magnitude  of  X^\-(02;c0|,‘0)j,(U2) 
for  each  liquid  was  measured  as  a  Ainction  of  probe  energy  for  a  range  of  wavelengths  between 
610  and  4S4nm  (2.04  to  2.73eV).  The  results  are  shown  in  Figure  8  where  the  scale  was 
determined  by  comparing  the  literature  value  fm  €$2  >  2. 14  x  10**^  esu)^^  to  the  measured 

ratio  of  the  optical  Kerr  constants  of  nitrobenzene  to  CS2  at  a  probe  wavelength  of  S14.5nm 
The  error  in  the  data  is  calculated  from  equation  (21)  by  attributing  a  5%  uncertainty  to  the 
centroids  €(<02)  ^  C(o\)  4nd  to  the  laser  powers  Fi„(oi>2)  &nd  Ppogp*  This  procedure  resulted 
in  a  overall  statistical  uncertainty  of  10%;  systematic  uncertainty  is  reduced  by  taking  repeated 
sets  of  dispersion  data.  Each  data  point  is  an  average  over  five  trials  that  thus  reduces  the 
uncertainty  to  5%. 

5.  DISCUSSION 

5.1  ELECTRONIC  EXCITATIONS 

As  completely  described  in  the  work  of  Orr  and  Ward,2®a  quantum  mechanical 
perturbation  treatment  of  the  electronic  third  order  susceptibility  results  in  an  analytic  expression 
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for  the  electronic  OKE  susceptibility  that  contains  terms  proportional  to  the  product  of  four 
transition  moments  and  inversely  proportional  to  a  product  of  three  energy  difference  terms. 
Near  resonance,  the  energy  denominator  of  a  single  term  can  dominate  the  frequency  dependence 
of  the  nonlinear  susceptibilty.^^  Under  these  condidons,  the  electronic  dispersion  of  the  third 
order  susceptibility  for  an  (o}|'KU2)  resonance  can  be  approximated  by  the  phenomenological 
expression 


OR(-(l>j;(0,,-ffl,,0)2)- - - - -  -  (30) 

-  0)^  -  <02  “ 

where  F  is  a  damping  constant,  the  sum  of  the  pump  and  probe  frequency,  (cO]-ko2)>  is  assumed 
Icxated  near  the  electronic  transition  frequency,  of  the  first  naolecular  excited  state,  and  all 
the  factors  that  are  frequency  independent  near  this  resonance  are  assumed  to  be  contained  in  the 
constant  A. 

The  measured  electroiuc  contribution,  au(-O)2;t0j,~<Di,(02)  as  a  function  of  probe 
frrequency  with  fixed  pump  frequency,  (0|,  then  follows  the  expression 


10^12 


A2 

((Oq.  C02)2+r2’ 


(31) 


where  -  (0|  «  A  major  advantage  of  the  pump-probe  experiment  is  that  the  measured 
resonance  does  not  truneh  the  frequency  of  either  the  pump  or  the  probe  beam.  When  the  sum  of 
the  pump  and  probe  frequency  matches  the  material  resonance,  no  single  photon  absorpnon 
occurs,  making  it  possible  to  probe  strong  resonances  with  little  depletion  of  either  beam,  and 
thereby  allowing  thorough  study  of  the  electronic  contribution  to  X^^\-<a2;c0i,“<0i.ti)2) 
centrosymmetric  molecules,  the  states  that  are  accessible  are  necessarily  two-photon  states  v*  hi«. 
have  been  shown  to  be  significant  contributon  to  the  hypcrpolarizability  of  orgaiuc  systems  - 
In  noncentrosymmetric  media  one-  and  two-photon  selection  rules  do  not  apply  and  any  state  s « 
priori  accessible.  When  (<0|+(a2)  is  far  from  any  material  excitation  energy,  the  contnbut »  « 
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from  all  of  the  nonresonant  electronic  excitations  will  be  refened  to  as  is  present  in  all 

materials  and  will  constitute  a  small  background  contribution  when  is  appreciable. 


5.2  MOLECULAR  REORIENTATION 

As  has  been  long  known,  molecular  reorientadon  can  occur  if  the  molecular  linear 
polarizability  is  anisotropic;  for  a  full  discussion  of  molecular  reorientation  the  reader  is  referred 
to  Wagner,  et  al.^  An  applied  electric  Held  can  cause  alignment  of  the  constituent  molecules, 
imparting  a  birefringence  to  the  ordinarily  isotropic  liquid.  The  energy,  W,  of  a  molecule  is 
given  by  dipole  coupling  to  the  field. 


mm  -  -P  E  -  -  a.xe,^)<E^E>, 

ij 


(32) 


where  E  is  the  pump  field,  P  the  induced  dipole  moment,  (Xjj  the  linear  polarizability,  and  the 
brackets  denote  a  time  average. 

If  we  assume  that  each  molecule  in  the  liquid  acts  independently,  then  for  uniaxial 
molecules  the  partition  function  can  be  written  as 
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exp  ^PAj(aj^+(aj-aj^)cos^)]  dcos8 , 
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where  and  a|  are  the  two  independent  components  of  the  molecular  linear  polarizability,  0  the 
Boltzman  factor  and  A^  the  magnitude  of  the  electric  field.  The  macroscopic  susceptibiliry  is 
given  by. 


a  ■  — 


i  g(Af). 


(34) 


The  refractive  index  can  be  expanded  in  a  Taylor  senes  and  the  first  two  terms  are 


-  ^  .  1  _or  -2 

n  *  Dfl  +  2  ’ 


(35) 


where  is  the  linear  refractive  index  and  n^*^  is  the  nonlinear  index  from  molecular  reorientation. 
The  expression  for  n^'  in  the  small  signal  approximation  for  the  optical  Kerr  effect  owing  to 
molecular  orientation  is^ 

where  N  is  the  number  density  and  T  is  the  temperature.  Cleary,  then,  n^'  will  adopt  the 
dispersion  of  the  polarizability  anisotropy  which  is  relatively  flat  far  from  one-photon 
resonances.  The  lattice  contribution  D  has  Fourier  components  at  zero  frequency  (molecular 
reorientation)  and  at  die  difference  frequency.  Am  (RKES),  and  is  conventionally  written  as 

D  =>  T  B(0)  +  { A(A<a)  +  j  B(Am) ) ,  (37a) 

and  we  can  write 

2nJ((02)nJ 

0(0)-—- - ■.  (37b) 

3jcFl 

where  is  the  local  field  factor,  and  B(0)  is  that  part  of  the  lattice  contribution  arising  from 

molecular  reorientation. 

The  classic  molecular  reorientation  effect  can  be  thought  of  as  the  d.c.  response  of  ihx 
molecular  nuclei  to  the  intense  applied  field;  there  is  also  the  possibility  of  resonances  u.ith 
vibrational  and  rotational  excitations  of  the  molecular  nuclei  expressed  as  a  resonant  e 
contribution  {A(Aci))  +5’B(A(a)).  When  the  difference  between  the  pump  and  probe  laser 
frequencies  is  near  a  Raman  resonance  a  RIK£‘^  is  observed,  quantified  by  the  Raman  hre  s 
frequency,  linewidth,  and  amplitude.  The  optical  Kerr  effect  can  be  used  to  measure  Ram  in 

'  ;6 


spectra  by  tuning  the  probe  fluency  coj  through  successive  resonances  or  by  using  a  broad 
band  probe  hence  getting  the  entire  spectrum  in  one  laser  shot  It  is  also  clear  that  by  measuring 
X^^\-o)2;cO|,-<a^,(a2)  near  a  Raman  resonance  one  can  measure  the  non-resonant  contribution  to 
through  a  knowledge  of  Raman  scattering  data.  In  this  way  accurate  values  of  X^^^  have 
been  obtained  for  various  Kerr  liquids;^^  a  complete  review  of  these  Raman  related  OK£ 
techniques  has  been  prepared  by  Levenson  and  Song.^  As  relates  to  the  current  studies,  it  is 
important  that  the  difference  frequency  002-  cOj  remain  far  from  fundamental  and  overtone 
Raman  active  modes  if  a  simple  interpretation  in  terms  of  sum*  ((D|  -f  axj)  and  one*  (coj.co^) 
photon  dispersion  is  to  be  possible.  Further,  a  general  approach  to  the  interference  between 
resonant  (X^^)  and  nonresonant  (X^)  terms  in  OKE  measurements  is  necessary. 

5.3  LOCAL  FIELD  EFFECTS 

The  macroscopic  susceptibility  is  related  to  the  microscopic  susceptibility  through 
frequency  dependent  local  field  factors.  While  there  are  various  prevalent  models  (Loren  tz 
Lorenz,  Onsager)  all  of  them  give  the  field  factors  in  terms  of  the  refiractive  index.  For  a  third 
order  process,  the  local  field  correction  facuv,  F^^,  is  a  produa  of  four  factors,  and  for  the  pump- 
probe  configuration  is  (rf  the  form, 

FL-fi(n(tDi))fi(n((02)).  (38) 

where  there  is  a  local  field  frntor  for  each  incoming  and  out-going  wave.  Since  only  two  of 
four  local  field  factors  are  frequency  dependent  in  this  pump-probe  experiment,  the  dispersion  ^ 
the  third  order  susceptibility  will  be  less  sensitive  to  the  local  field  facton  than  expenmc'-s 
where  all  frequencies  are  varied  as  in  third  harmonic  disperison  measurements  or  in  two  bei.- 
experiments  where  pump  and  probe  beams  are  derived  from  the  same  source. 

In  typical  dispersion  experiments,  the  range  of  the  probe  wavelength  is  about  200  -  -- 
The  change  in  index  over  this  range  is  typically  less  than  5%  for  our  samples  and  this  resu  :  s 
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changes  of  from  Lorentz-Lorenz  field  factors  (f^  “  ^(n^+2)(n^+2))  by  about  9%.  This  is 
small  compared  to  the  effects  of  dispersion  on  resonance  that  typically  enhance  the  signal  by  a 
factor  of  two  or  nx»e  in  states  of  weak  oscillator  strength  that  are  at  the  detection  limit  of  linear 
absorption  measurements.  For  strong  resonances,  the  dispersion  from  the  local  field  factors  is 
negligible.  For  the  remainder  of  this  paper  we  will  assume  that  the  local  field  variation  in  the  the 
nonlinear  susceptibility  is  negligible  over  the  frequency  range  interrogated,  so  the  local  field 
factor  will  not  be  explicitly  shown.  However  if  accurate  comparisons  of  different  materials  is 
required  the  local  field  factors  must  be  included. 

4 

5.4  INTERFERENCE  BETWEEN  RESONANT  AND  NONRESONANT  PARTS  OF 

The  pump  and  probe  laser  frequencies  used  in  the  present  measurements  are  far  off 
resonance  for  both  benzene  and  nitrobenzene  with  little  dispersion  of  the  linear  refractive  index 
(~2%)  over  the  same  range  (see  table  1).  Ifence,  any  dispersion  in  the  measured  signal  is  due  to 
electronic  sum-  and  difference-resonances  as  evident,  for  example,  in  Equation  (31).  If  the  sum 
frequency  corresponds  to  a  weak  electronic  excitation,  the  total  susceptibility  will  be  a  sum  of 
resonant,  Xr  \  and  nonresoiumt,  X^,  components, 

+  XgJ.  (39) 

where  the  magnitude  is  given  by 

IX^^^  -  [  IXg¥  +  +  2RE[Xg’ 

and  the  spectrum  will  exhibit  the  interference  between  the  resonant  and  noiuesonant  portions  xs 
first  studied  extensively  by  Levenson,  et  aL“  In  terms  of  our  previously  defined  parameters 


-  I  <»R  +:^7l{2A„(A(0)  +  B„(Aa»). 
n  n 
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(41b) 
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where  is  the  resonant  contribution  to  the  hyperpolarizability  centered  at  frequency  Q^g,  and 
A„(Aci))  and  B„(Aco)  are  the  resonant  Raman  contributions  to  the  nonlinear  susceptibility  at  center 
frequency  Q)„.  Thus  knowledge  of  can  allow  determination  of  the  separate  resonant 
contributions  while  a  prediction  for  X^^say,  from  molecular  orbital  calculations)  can  be  used  to 
determine  the  nonresonant  nonlinear  susceptibility. 

The  nonresonant  susceptibility  X||^  has  a  weak  frequency  dependence  that  can  be 
observed  if  measurements  are  marfe  over  a  broad  enough  frequency  range.  and  B(0)  must 
be  treated  separately,  since  the  dispersion  of  B(0)  follows  that  of  the  linear  polarizability 
anisotropy,  while  that  of  is  related  to  the  second  hyperpolarizability.  One  immediate 
ramification  is  tiiat  B(0)  can  only  derive  frequency  dependence  from  a  product  cf  two  one*photon 
energy  denominaton,  while  0^  obtains  dispersion  from  a  product  of  three  one*  and  two-photon 
energy  denominaton. 

Several  methods  exist  for  developing  such  an  expansion  for  ranging  from  a  direct 
expansion  of  the  time  dependent  perturbation  theory  result  to  calctilations  based  on  isotropic 
anhannonic  oscillator  models  of  the  nonlinear  susceptibility.^  In  the  model  of  Owyoung,^^  Ovr 
is  written  as 


<yg5^(-<02;<i02,o»i,-<»i)  -  a^ca2)a2(©i),  (42) 

where  a  is  the  isotropically  averaged  linear  polarizability.  If  we  suppose  we  can  expand  .‘le 
linear  polarizability  according  to  its  representation  as  a  harmonic  oscillator 


a(0)_)  » cXo(l4-((i)^Qo)2  + 


(43) 


where  is  the  firequency  of  a  large  oscillator  strength  high  energy  excited  state,  then  we  get  the 
following  result  for 

(1-K^)^)  (l+(^)^)^.  (44) 

If  we  let  ctj  (1  +  =  Oq  then 

Onr  =  <Jo(1+ 2(^)2),  (45) 

so  the  (O2  dependence  is  quadratic  to  lowest  order. 

Alternatively,  we  can  use  the  result  for  a  from  third  order  perturbadon  theory  in  the 
dipole  approximation.  If  damping  is  assumed  to  be  negligiUe  then  it  can  be  shown  that  this  must 
be  an  expansion  in  even  powers  of  the  optical  frequency.  If  damping  is  present  as  when  our 
probe  frequency  approaches  a  sum-resonant  frequency  then  the  expansim  will  contain  terms 
linear  in  0)2*  Explicit  knowledge  of  molecular  wavefunctions  and  energy  levels  would  be 
necessary  to  obtain  the  coefficients  in  this  expansion  so  we  simply  write 

a  Cj  +  C2®®2’ 

From  (37),  B(0)  obtains  dispersion  from  the  polarizability  anisotropy.  Since  a  single  one 
photon  energy  denominator  contributes  in  this  case  we  might  expect  B(0)  to  have  less  dispersion 
than  however,  B(0)  is  generally  larger  than  in  liquids,  so  the  frequency  dependence  't 
B(0)may  be  appreciable.  For  the  sake  of  simplicity  we  assume  B(0)to  be  frequent  v 
independent  within  the  experimentally  nreasured  regions,  and  discuss  the  success  of  this  ansaa  i 
posteriori. 

In  terms  of  equations  (41)  we  can  write 


130 


(47) 


-  5»«»NR  +  +  B„(Acd)}+  I  . 

n  n 

Whichever  model  we  adopt  there  will  be  certain  shared  characteristics  in  the  spectra,  notably  the 
appearance  of  Fano-like^  resonance  peaks  when  0)^'K02  =  ,  resulting  firom  the  interference 

of  X^^^and  X When  X^^  is  a  Lorentzian  resonance,  the  upper  portion  of  Figure  9  shows  the 
absolute  magnitude  of  the  susceptibility  with  X^  a  0.  In  the  bottom  portion,  X^j^  is  given  by 
either  (48)  or  (49),  which  are  valid  if  other  electronic  excitations  are  far  from  resonance. 

5.5  X^^^  OF  BENZENE  AND  NITROBENZENE 

The  observed  dispersion  for  nitrobenzene  can  be  understood  in  terms  of  the  gas  phase 
linear  absorption  spectrum  which  gives  the  energies  and  oscillator  strengths  of  the  electronic 
excited  states  in  the  near  ultraviolet^*^  The  general  rise  in  the  |  X^^  |  data  (Figure  8)  may  be 
related  to  the  onset  of  a  sum  resonance  with  the  high  oscillator  strength  state  at  5.11  eV  (gas 
phase),  while  the  fine  structure  at  a  probe  energy  of  2.6  eV  is  associated  with  a  weak  gas  phase 
transition  at  4.35  eV  .  In  neat  solution  the  5.11  eV  band  is  red  shifted  by  0.65  eV,^’  thus 
shifting  the  sum  resonances  by  an  equivalent  amount  The  rather  complex  fine  structure  may 
result  from  interference  between  a  sum  resonance  with  the  red-shifted  4.35  eV  band  and  a 
difference  resonance  (<02  ~  ^1)  ^  infrared  overtone  band  centered  at  1.48  eV.  The 

data  was  analyzed  using  (47)  with  Lorentzian  resonances  centered  at  probe  energies  of  2.S5  e  v 
(sum  -  electronic)  and  2.66  eV  (difference  •  nuclear)  and  experimentally  determined  widths  ot 
0.09  eV  and  0.04  eV,  respectively.  If  was  assumed  to  have  no  frequency  dependence,  a 
good  fit  could  be  achieved  to  all  but  the  point  at  the  highest  probe  energy  (Figure  8).  The  2  55 
eV  resmiance  amplitude  was  15%  of  the  nonresonam  background,  while  the  2.66eV  peak  ^  as 
weaker  with  an  amplitude  of  9%  of  the  background  and  the  estimated  error  in  these  amplitudes 
was  approximately  5%.  Far  from  resonance  this  model  predicts  £  0*5  constitutes  15%  ot  -e 

m 

nonresonant  X^\  and  this  is  in  good  agreement  with  other  experiments.^^ 


If  allowed  to  have  fluency  dependence,  then  both  models  (45)  and  (46)  predict 

that  B(0)  must  be  negative,  and  this  is  ruled  out  by  the  posidve  definite  expression  (36).  Since 
the  next  strong  sum  resonant  excited  state  is  at  a  probe  energy  of  3.3  eV,  it  may  be  close  enough 
to  prevent  the  binomial  expansions  (45)  and  (46)  from  being  applicable,  thus  suggesting  that  a 
three  Lxirentzian  model  might  be  necessary ;  this  was  not  attempted  because  there  was  no  data  at 
higher  energies.  Furthermore,  the  large  resonant  contribution  may  completely  mask  the 
nonresonant  frequency  dependence,  including  any  dependence  from  B(0)  that  we  have  neglected 
A  likely  reason  for  the  failure  of  the  expansion  (45)  to  explain  the  nonresonant  behavior  is  that 
the  coefficient  of  the  quadratic  term  as  predicted  by  the  theory  is  incorrect,  as  has  been  observed 
in  dc-field  induced  second  harmonic  generation  experiments.^^ 

Benzene  possesses  primary  excited  states  at  4.75, 6.14  and  6.74  eV  in  the  gas  phase,  and 
these  suffer  from  little  dielectric  shift  in  the  liquid.  The  I I  spectrum  of  benzene  (Figure  8) 
shows  less  pronounced  structure  than  that  of  nitrobenzene  because  of  the  smaller  dielectric  shift 
and  the  preservation  of  the  parity  selection  rule  in  the  centrosymmetric  benzene  molecule.  The 

data  is  largely  dispersion-free  aside  from  some  structure  at  a  probe  energy  of  ~  2.5  eV 
that  coincides  with  a  weak  optical  transition  near  3.7  eV  reported  for  liquid  benzene. The  first 
strong  excited  state  at  4.86  eV  is  not  experimentally  accessible  because  of  the  small  solvent  shift 
and  the  lack  of  a  sum-photon  contribution  from  that  state.  We  note  that  the  measured 
nonresonant  value  of  I X^^  i  benzene  at  Xj  >  514.5  nm  was  about  two  times  less  than  that  of 
nitrobenzene,  which  is  in  agreement  with  the  literanire  when  laser  pulse  length  corrections  are 
accounted  for.^^ 

As  in  the  case  of  nitrobenzene,  the  benzene  spectrum  shows  the  structure  that  results 
from  the  interference  between  a  Lorentzian  resonance  and  the  background  of  molecular 
reorientation  and  non-resonant  electronic  transitions.  If  the  polynomial  expansion  (^6)  s 
employed  the  benzene  spectrum  can  be  well  desenbed  by  a  single  Lorentzia  ■  with  A=-0  Ob’S  < 
10'‘^  esu-eV,  ti)„  »  2.50  eV,  T  =  0.01 1  eV.  C,  +  B(0)  =  24.9  x  10-‘3  esu  and  q  -  -4.28  x  1 0 
esu/eV  with  estimated  errors  as  above.  The  Lorentzian  parameters  are  reasonable  for  a  '''■ei* 


oscillator  strength  transition  in  a  room  temperature  liquid  and  describe  the  electronic  contribution 
of  the  weak  state  at  3.7  eV  in  benzene;  furthermore  they  suggest  that  this  state  has  appreciable 
sum-resonant  character.  The  sign  of  is  unexpected  and  may  result  from  frequency 
dependence  in  B(0)  or  from  angular  cotreladons  between  adjacent  benzene  molecules,  which 
have  not  been  considered.  The  total  electronic  contribution  can  only  be  determined  if  the  probe 
frequency  is  tuned  through  the  significant  sum-resonant  sutes  and  the  Lorentzian  parameters 
determined  as  above.  For  the  present  study,  we  illustrate  the  technique  by  determining  the 
electronic  contribution  from  one  weak  excited  state,  while  a  and  B(0)  remain  unseparated. 

6.  CONCLUSION 

An  experimental  and  theoredcal  study  of  OKE  dispersion  in  isotropic  media  has  been 
presented.  A  theoredcal  analysis  based  on  Maxwell's  equadons  provides  the  connection  between 
experimentally  measured  quantities  and  Xjjy(-(i)2;o)j,-<i)|,CD2)«  An  optimized  pump-probe 
depolarization  technique  was  developed,  and  conventional  theoretical  approximations  carefully 
studied.  Thresholds  for  pump  saturation  and  thermal  leasing  were  obtained,  and  experimental 
powers  accordingly  adjusted.  The  OKE  dispersion  data  is  used  to  determine  the  resonant  and 
nonresonant  contributions  to  ?^^(-<Q2;o>i,-0)|.n>2)  near  sum-frequency  resonances,  and  can  be 
used  to  separate  the  nonresonant  electronic  and  lattice  contributions  when  a  broad  enough 
frequency  range  is  investigated.  The  data  were  analyzed  according  to  a  simple  model  that 
required  comparison  with  Hne***  absorption  spectra,  a  polynomial  expansion  of  the  nonresonant 
electronic  contribution,  and  consideration  of  dipole  selection  rules.  The  complex  nitrobenzene 
spectrum  was  seen  to  result  from  an  interference  between  a  sum-frequency  resonance  (to,  >  i 
a  difference- frequency  resonance  (0)2  ~  Ot),  and  nonresonant  contributions.  A  simple  model 
succeeded  in  explaining  the  results  for  benzene  for  which  a  single  weak  electronic  transiaon  it 

3.7  eV  was  probed  far  from  any  other  resonances  in  the  presence  of  a  nonresonant  background 

\ 

The  general  formalism  presented  in  this  study  can  be  adopted  for  the  study  of  resonant  X  ^ 
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resonant  processes  in  any  other  media  of  interest,  and  the  experinsental  data  can  be 
compared  to  a  microscopic  theoy  based  upon  many  electron  molecular  orbital  calculations.^ 


nCURE  CAPTIONS 


Figure  1.  Experimental  layout  where  HWPl  is  a  half  wave  plate  at  1.06iim,  PI  and  P2 
polarizers,  A1  and  A2  analyzers,  BSl  and  BS2  beamsplitters,  MDMl  and  MDM2  motor  driven 
mirrors,  KDP  mixing  crystals,  PD  photodiodes,  FI  and  F2  filters  and  PMTs  arc  photomultiplier 
tubes. 

Figure  2.  Typical  pump  and  probe  beam  profiles.  The  pump  profile  is  seen  to  agree  with  the 
parabolic  approximation  (dashed  line)  within  the  probe  beam  waist 
Figure  3.  Poissonian  pulse  height  distribution  of  normalized  PMT  signal  accumulated  over 
30,000  laser  shots. 

Figure  4.  Log  plot  of  the  heating  decay  curve  of  nitrobenzene  with  a  Hirie  constant  of  1.42  sec. 

Figure  5.  Thermal  defocusing  profile  in  nitrobenzene  as  determined  with  a  «^*nning  pinhole. 

Figure  6.  Argon  ion  beam  profiles  at  a  wavelength  of  488  nm  and  a  power  of  60  mW  through 
nitrobenzene  (solid  curve)  and  without  sample  (dashed  curve). 

Figure  7.  Kerr  signal  as  a  function  of  half  wave-plate  angle  (see  Figure  1)  in  the  pump  beam 
path.  Curve  represents  the  expected  sin^28  dependence. 

Figure  8.  I |  dispersion  data  for  nitrobenzene  is  given  by  the  solid  points  and  benzene  data 
is  given  by  the  open  points. 

Figure  9.  The  upper  curve  shows  a  Lotentzian  form  of  the  absolute  value  of  the  susceptibtbty 
for  a  two  level  modeL  while  the  bottom  curve  shows  the  behavior  of  the  absolute  value  of  an 
offset  Lotentzian. 
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Figure  2.  Typical  pump  and  probe  beam  profiles.  The  pump  profile  is  seen  to  agree  with  the 
parabolic  approximation  (dashed  line)  within  the  probe  beam  waist 
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Figure  6.  Argon  ion  beam  profiles  at  a  wavelength  of  4«8  nm  and  a  power  of  60  mW  through  | 

nitrobenzene  (solid  curve)  and  without  sample  (dashed  curve). 
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I^Surt  7.  Ken  signil  as  a  function  of  half  wave*plate  angle  (see  Figure  1)  in  the  pump  >vam 
path.  Curve  represents  the  expected  sin^29  dependence. 
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Figure  8.  1 I  dispenion  data  for  nitrobenzene  is  given  by  the  solid  points  and  benzene  data 
is  given  by  the  open  points. 
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